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I.  INTRODUCTION 


The  determination  of  the  two-point  velocity  correlation  tensor  R^j  is  a 
major  aim  of  our  research  because  R^  contains  full  information  on  the 
Reynolds  stress  tensor,  the  scale  andJthe  structure  of  an  arbitrary  turbu 
lent  flow.  The  tensor  Rjj  Is  defined  by 

Rij(x,y)  =  <(uj(x)uj(y)^ 


We  propose  to  determine  R-jj  by  solving  a  closed  equation  for  this 
quantity  (Refs.  1,  2).  Closure  requires  modeling  of  the  nonlinear  velocity 
correlation  terms  analogous  to  those  now  familiar  in  second-order  closures 
for  the  Reynolds  stresses  (Ref.  3).  Because  Rij  contains  scale  and 
structure  Information,  closure  of  the  Rij  equation  obviates  the  need  for 
Introducing  ad  hoc  information  on  these  features  of  the  turbulence. 

As  a  first  attack  on  the  integration  of  the  modeled  Rij  equations  for  a 
general  flow,  we  simplify  the  integration  by  introducing  the  approximation 

*u-#V|,|)  (1) 


which  is  the^fir^t  t£rm  in  the  moment  expansion  of  Rii  in  the  relative 
coordinate  r  *  x  -  y  .  Here  Aij  is  the  angular  or  Batchelor  average  of 
the  integral  scale  for  the  appropriate  component  of  the  two-point  correla¬ 
tion  tensor;  i .e. , 


(2) 


This  "first  moment"  approximation  yields  coupled  equations  for  the  Reynolds 
stress  and  for  the  scale  tensor.  Since  both  quantities  are  single-point 
variables,  the  approximation  Is  relatively  simple  when  compared  to  the  full 
system,  which  is  a  two-point,  multicomponent  equation. 


We  note  that  for  Isotropic  turbulence,  as  discussed  in  Ref.  1,  we  have, 
with  A  =  A|ck/3 


1.  Sandri,  G.,  "A  New  Approach  to  the  Development  of  Scale  Equations  for 
Turbulent  Flows,"  A.R.A.P.  Report  No.  302,  April  1977. 

2.  Sandri,  G. ,  "Recent  Results  Obtained  in  the  Modeling  of  Turbulent  Flows 
by  Second-Order  Closure,"  AF0SR-TR-78-0680,  February  1978. 

3.  Donaldson,  C.duP.,  "Construction  of  a  Dynamic  Model  of  the  Production 
of  Atmospheric  Turbulence  and  the  Dispersal  of  Atmospheric  Pollutants," 
Workshop  on  Micrometeorology  (D.A.  Haugen,  ed.),  American  Meteorologi¬ 
cal  Society,  Boston  (1973)7 PP»  313-392. 


A 


f(r)  dr 


(3) 


f  A  =  \  \k  “  Lf  *  f 

Jo 

where  r  =  |r|  .  Since  Lf  is  6.1.  Taylor's  integral  scale,  Eq.  (3) 
justifies  our  statement  that  the  tensor  Aij  represents  an  integral  scale. 

The  first  moment  approximation  is  based  on  the  assumption  that  the  size  of 
the  eddies  is  small  compared  to  the  scale  of  the  mean  quantities.  This 
assumption  is  not  usually  valid,  so  that  for  a  number  of  flows  we  must 
consider  the  full  R-jj  equation.  Fortunately,  however,  we  can  test  the 
first  moment  approximation  quite  stringently.  In  fact,  (i)  the  first  moment 
approximation  is  valid  for  homogeneous  flows,  whether  sheared  or  not,  because 
for  such  flows  the  scale  of  the  mean  is  infinite,  and  (ii)  several  experi¬ 
ments  have  been  performed  on  homogeneous  flows,  with  and  without  shear, 
which  yield  information  on  several  turbulence  quantities. 

Our  model,  which  contains  only  two  adjustable  parameters,  is  required  to 
reproduce  the  main  features  of  both  grid  and  sheared  turbulence  data.  We 
find  a  satisfactory  overall  fit  with  parameter  choices  that  are  used  in 
second-order  closure  for  general  flows.  The  calculated  integral  scales  for 
sheared  turbulence  are  strongly  directional.  Thus  the  model  allows  one  to 
compute,  for  the  first  time,  a  first  approximation  of  the  eddy  structure  in 
a  sheared  turbulent  flow. 


Proceeding  a  step  beyong  the  first  moment  approximation,  we  show  that  the 
familiar  dissipation  equation 


-2  n3 

m — 2b 


can  be  deri ved  from  a  simple  model  of  the  "cascading"  term  (the  divergence 
of  the  triple  velocity  correlation) . 
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II.  DERIVATION  OF  CLOSED  EQUATIONS  FOR  THE  TWO-POINT  CORRELATION  TENSOR 


We  start  with  the  Navier-Stokes  equations  for  the  velocity  field  u-j  and 
the  kinematic  pressure  p  (pressure/density) 


3u<  ** 

uk  3^  +  3x7 


2 

v7£ui 


(4) 


0 


(5) 


The  velocity  and  pressure  are  decomposed, foil  owing  Reynolds,  into  mean  and 
fluctuating  parts 


ui 


ui 


(6) 


P  ■  P  +  P' 


(7) 


Substituting  (6) 
for  the  mean  and 
equations  can  be 

RU 


and  (7)  into  (4)  and  (5),  one  obtains  separate  equations 
fluctuations.  After  multiplication  and  averaging,  these 
cast  into  the  form  of  equations  for  the  correlation  tensor 


R-jjCx.y)  *  <(uj(x)uj(y£> 


(8) 


namely. 


(9) 


3 


and 


(10) 

(ID 


It  is  convenient  to  recast  these  equations  by  introducing  "centroid"  and 
"relative"  variables  as  follows: 


x 

it 


(12) 


(13) 


Using  the  chain  rule  for  differentiation,  as  detailed  in  Ref.  1,  the  dynami¬ 
cal  equations  for  become 


3R, 


at 


uk(x)  +  uk(y)  a 


T 


3  u .  3u . 


+  (“k^>  -  'ukl*>% 


Rij  + 


R.  . 


Rlk  d*ck  3xck 

-  [<(^(x)u£(x)uj(y))>  +  <^u:(x)u^(y)u^(y)^>] 

ck 


—  [^u!(x)u£(x)uj(y)^  -  <^C*)u£tf)ultf)>] 


|^<pTx)u^)>  +  |^<uj(*)Pltf)> 


cj 


/  ,-k  3ui^)\  /3u)(x)  \  ' 

Kp(x,4r) +  <^rp,(5,)j 

2  /3u{(x)  3uUy)\ 

^’V2',v4r^> 


(14) 
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and  the  continuity  equations  become 


(15) 


!!i I 
9xcj 


(16) 


where  the  arguments  in  the  correlations  are  given  by  the  inverse  of  Eqs.  (12) 
and  (13),  namely. 


x 


(17) 


y  *  xc  +  £  (is) 

For  homogeneous  turbulence,  the  derivatives  with  respect  to  the  centroid 
vanish  when  acting  on  any  correlation  but  not,  in  general,  on  the  mean 
velocity  or  on  a  fluctuation.  We  then  have,  for  homogeneous  turbulence,  the 
dynamical  equations 


3Uj  3u.j 
R,k  ^ 4  ^  R« 


+ 


§f-  <u!(x)u,[(x)uj(y))>  -  <(uj(x)u£(y)uj(y)]> 


’  _  3u ;(y) 

P’(X)  -r-J- 


3u!(x)  „ 

1 _  « I  IT. 


3xci 


3x 


p'(y) 


CJ 


/3ui(x)  3u-(y)' 


Sxck  axck 


and  the  continuity  conditions 
3^i  j 

— *  0 
3ri 


(19) 


(20) 

(21) 


The  terms  in  Eq.  (19)  that  prevent  closure  of  the  equation  are  the  last 
three.  In  order  to  achieve  closure,  we  Introduce  generalized  transport 
models  as  follows: 


(i)  Dissipation  equation.  We  shall  prove  in  Section  V  that  a  model  for 
dissipation  "is  not  needed  once  a  model  for  eddy  size  rearrangement  is 
specified.  A  full  specification  of  this  latter  is  not  available  at  present; 
therefore,  we  provisionally  maintain  a  model  for  dissipation  which  is  given 
by 


-2v 


/  3u!(x)  3ul(y) \ 
\^dT  3xck  / 


(22) 


( i i )  Intercomponent  rearrangement.  We  give  a  general ization  of  the  famil¬ 
iar  tendency- towards- isotropy  model  by  introducing  a  tensor  Cy  that 
guarantees  that  continuity  is  satisfied. 


3Uj(y)  ^  3ul(x) 


3x 


Cl 


3x 


CJ 


R,,  -  4  .R  + 

lj  3  l]  aa 


(23) 


The  tensor  C-j j  is  assumed  to  satisfy 


Ci(j(r  -  0)  -  0 


(24) 


(25) 


These  properties  are  utilized  in  Section  III  to  derive  the  Reynolds  stress 
and  scale  tensor  equations.  In  the  Appendix,  we  prove  a  rather  remarkable 
property  of  C-j j  for  isotropic  turbulence.  If  Sy  satisfies  continuity 
(as  we  assume  here),  the  tensor  Cy  is  uniquely  determined  by  continuity 
and  it  is  such  as  to  make  the  pressure-velocity  correlations  vanish  as  is 
necessary  for  isotropic  turbulence.  A  proof  of  the  latter  statement  is 
given  in  Ref.  4,  p.  51. 


(iii)  Eddy  size  rearrangement.  In  addition  to  the  two  generalizations 
given  above,  we  need  a  model  for  the  spatially  homogeneous  part  of  the 
triple  velocity  correlation  which  represents  the  nonlinear  effects  of  local 
turbulent  convection.  These  nonlinear  effects  correspond  to  either  eddy 
break-up  (cascading,  when  wave  vectors  add)  or  merging  of  eddies  (when  the 
wave  vectors  subtract).  We  assume  the  model 


lr-  -  <<i|(x)U|;(y)Uj(y)>]  =  V  a  [Rjj  .  S)jj  (26) 


4.  Batchelor,  6.K.,  The  Theory  of  Homogeneous  Turbulence,  Cambridge 
University  Press,  I55TT 
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where  S. .  satisfies 

1  J 


(27) 


S 


ij 


(r  *  0) 


d  r  =  0 


(28) 


4irr 


A  simple  example  of  Ni  ■  is 

■  J 


where  n  is  a  positive  integer.  Substituting  the  models  (22),  (23)  and 
(26)  into  (19)  yields 


(29) 


with  (20)  and  (21),  i.e.,  continuity,  holding. 


We  also  introduce  the  spectral  tensor  $jj  in  order  to  present  the 
equations  in  separation  space  as  well  as  wave  number  space.  The  definition 
is 

♦ijtf)  *  — f  ^Ar)  e1***  dr 

(2tt)3  J  1J 

The  unmodeled  equation  for  the  spectral  tensor  is 


7 


■  U'(*125jl  *  S21*2j)  "  1 


1  [k1PJ<^>  -  kJP1<-k>]  ’ 


where  we  have  Introduced  the  two  spectra 


pj(k)  ■  J*<^p,(x)uj(y)>  e‘ik‘r  ^L. 


and  the  triple  velocity  correlation  tensor 


Tiotj(r)  *  <u1(x)uo(x)uj(x  +  r))> 


Continuity  requires  the  Poisson  equation 


2  Mk)  -  M-M  r  -I 

tk  }  2  j - 2U'kl»2j  +  1ka[‘«<k>**aj‘-k)J 

Anticipating  the  development  of  Section  V,  we  note  that  the  following 
simple  model  gives  a  tolerable  picture  for  the  energy  spectrum  of  grid 
turbulence 


-1  » 

*  a  +  B  k  If  *i j *  2*ij  4  c  (k2  “?  ♦u  +  4k  w  *i j +  2* 

with 

A  *  0.075  ,  B  =*  0.615  ,  C  *  0.27 
which  have  been  adjusted  to  give  exactly  a  5/3  law  for  equilibrium. 


A  full  matching  of  the  energy  spectrum  as  a  function  of  time  is  not 
achieved  with  this  model;  mainly  the  peak  value  is  followed  correctly.  We 
attribute  this  feature  in  part  to  the  fact  that  a  grid  flow  Is  not  simply 
an  isotropic  flow.  In  Section  V  we  discuss  a  simpler  form  of  the  model  to 
demonstrate  an  important  relation  between  dissipation  and  eddy  size 
rearrangement. 
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III.  EQUATIONS  FOR  THE  REYNOLDS  AND  SCALE  TENSORS 


We  now  exploit  the  definition 


in  order  to  obtain  the  rate  equation  for  the  Reynolds  stress  (divided  by 
p  ).  We  let  r  -*>  0  in  (29)  and  find,  using  (24)  and  (27) 


It  »i“j  ■ 


- du,  au,  - 


-*(¥)W 


(30) 


The  scale  tensor,  A^j  ,  is  defined  in  terms  of  R^j  by 


(31) 


To  obtain  a  rate  equation  for  the  scale  tensor  Afj  that  includes  the  first 
order  information  of  turbulence  structure,  we  expand  Rij(r)  in  terms  of 
its  moments  In  pure  relative  coordinates  and  retain  the  lowest  term.  We 
obtain 


R1j(0 


V  A1j6d?l> 


(32) 


where  the  Dirac  function  of  the  magnitude  of  f  is  related  to  6(?)  by 

4(f)  .  Silfli  (33) 

4er* 

We  use  the  notation  r  -  |r|  .  The  approximation  (33),  as  we  remarked  in 
the  Introduction,  Is  best  justified  for  homogeneous  turbulence  because  In 
this  case  the  spatial  scale  of  the  mean  flow  Is  Infinite. 


We  now  apply  to  (29)  the  operation  f  dr/4irr*  and,  using  (25)  and  (28),  we 
find  J 


r 

! 


*  ?  [Aj  •  i  6ijq2Akk]  •  2tK*36u 


(34) 


Convection  of  A<h  does  not  occur  for  homogeneous  turbulence.  To  see  that 
this  fact  Is  a  consequence  of  the  moment  approximation  (32),  we  note  that 
(33)  gives 


where  we  performed  an  Integration  by  parts  and  used  continuity.  We  obtain 
the  final  form  of  the  equation  for  the  scale  tensor  by  substituting  Into 
(34)  the  formula 


3_ 

at 


*  at 


a_2 

Au  If- 


(35) 


and  the  contraction  of  Eq.  (30) 


i<)2  .  ii  8Ul  3 

H — z  “i“i  3*^7  - V 


(36) 


This  latter  Is  the  energy  equation.  The  result  Is 


12 


> 


Note  that  Eq.  (37)  us  Ibtalned  by  dividing  by  q£  and,  therefore.  It  Is 
meaningless  In  the  absence  of  turbulence. 

To  close  the  pair  (30)  and  (37),  we  choose 

A  *  5  \k  < 


The  reason  for  this  choice  was  discussed  In  a  previous  report  (Ref.  1). 
Contraction  of  (37)  and  use  of  (38)  yields  the  equation  for  the  mean  scale 


3A  .  2  . 

It  IA1k 


2ujui  3u. 

♦  v'q 


It  Is  Interesting  to  note  that.  In  this  model,  the  coefficient  of  the 
"production"  term  Is  not  a  universal  constant  and  receives  generally  compet 
Ing  contributions  from  the  Reynolds  stress  and  from  the  tensor  scale. 


IV.  SOLUTIONS  OF  THE  COUPLED  EQUATIONS 
FOR  STRESS  AND  SCALE  TENSORS 


In  this  section  we  give  two  analytic  solutions  to  the  coupled  equations  for 
the  stress  tensor  and  scale  tensor  equations.  We  use  subsections  to 
separate  the  different  calculations. 

A.  Equations  In  Standard  Coordinates 

The  centroid  vector  and  mean  velocity  vector  are  taken  to  have  components 


(x,y,z)  ,  (U{y) ,0,0) 

with  3U/3y  *  U'  *  constant.  The  relevant  components  of  the  Reynolds  stress 
equations  are  obtained  from  Eqs.  (30)  and  (36).  We  drop  primes  on  the 
fluctuations  and  give  a  form  useful  for  numerical  Integration  In  which  u? 
and  Ajj  are  calculated  from 


7 


2  ~l 


hi  *  3A  "  *22  "  A33 


(40) 

(41) 


The  other  relevant  components  of  the  stress  and  energy  equations  are 


It  “2  ’  I  (1  '  2b)  ak”  1  “2 


(42) 


If*?  * V- 1*? 


(43) 


If  u1u2  *  *u2  U'  ’  I  ulu2 


(44) 


u.u,  „2 


(45) 


For  the  tensor  scale  components,  we  obtain,  using  Eqs.  (37)  and  (39) 


ft  *22  u  ‘  1*22  +  (1  '  2b)q 


(46) 
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(47) 


ft  A33  "  '  T  A33  +  {1  ‘  2b)q 

Tt  A12  T  A12  ■  UA22 

a  uiu2  2 

|FA.2aj2U'A-|A12U’  +  v', 

where 

i  u.u,  _ 

-2  U'  ♦  J  (1  -  2b  -  v') 


(48) 

(49) 

(50) 


B.  Solution  of  the  Shearless  Equations 

Setting  U'  *  0  ,  we  see  that  equations  for  q  and  A  decouple  from  the 
tensor  components.  Introducing  the  deviators 


and  the  time 


We  have  the  set 


Integrating  (56)  as 


■  Vj  -  j  V2 

(51) 

°1j 

*  Au  •  i  6ij\k 

(52) 

X 

*  A 
q 

(53) 

ft  q“ 

b  3  A  «  V'  A 

■  x  q  •  at  A  x  A 

(54) 

ft  d1J  * 

b  ,  in  _  1  -  2b  -  v'  „ 

x  d1j  *  3t  °1j  x  °ij 

(55) 

ft T  “ 

b  +  v' 

(56) 

as 

T  ■ 

(b  +  v')(t  -  tj  ♦ 

^  A 

(57) 
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we  see  that  q  ,  A  ,  d^j  and  are  suitable  powers  of  (q  /T  )t  ;  for 
example,  0 


9  *  9o 


d1j  *  V 


A  «  A. 


(b  +  V)j&(t-tJ*l 
0 


(b  +  V')  /  (t  -  t)  +  1 
0 


T-b/(b+v‘) 

1/ (b+v  * ) 

-i  -v'/(b+v' ) 


(b  +  v 1 )  r  {t  ■  to)  +  1 


(58) 

(59) 

(60) 


A  good  fit  to  experimental  data  on  the  decay  of  turbulent  energy  and 
the  growth  of  eddy  size  for  grid  turbulence  Is  obtained  if  one  chooses 

t.  :  i  .  v'  :  0.075 


(61) 


We  then  see  that  for  large  times 

2  2 
"i**  ^ 


%  [•  2)f  ‘ 


'  du 


,[■ 


2q°t 

2  rt 


-5/4 

(62) 

r 

(63) 

which  shows  that  the  deviator  decays  with  a  power  about  four  times  larger 
than  the  energy.  This  substantial  difference  may  eventually  be  checked  In 
our  anisotropic  grid  flow. 

From  the  solutions  given  above,  we  can  verify  that  statistics  are  preserved 
by  the  model  equations  If  the  model  parameters  satisfy  certain  bounds.  We 
first  show  that  the  two  tensors  ujuj  and  A^j  are  positive  definite  from 
their  definitions.  Consider  an  arbitrary  (constant)  Ai  then. 


A1u1uj*j  *  (U‘A)^  1  0  (®4) 

the  equality  sign  holding  for  A  £  0  only.  Thus  uTuT  Is  a  positive 
definite  tensor.  1  J 

From  the  definition  of  the  scale  tensor  given  by  Eq.  (31),  using  Fourier 
transform  on  , 
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_2 

VAij 


where  the  power  spectrum  tensor  <Hj  Is  positive  definite  by  Khlutchlne's 
theorem.  Thus , 


VAiAijAj  i 


Thus,  Ai i  Is  positive  definite  because  q*  Is  positive  as  a  consequence 
of  (64).  J 

Using  the  solution  (61)  and  an  analogous  solution  for  Ajj  ,  we  find 

/„  \  -l/(b+v' ) 


U^U)  *  Ui  Uj  (°)  ^ 


lx  „ 2  % 


-2b/(b+v' ) 


ft-) 


-l/(b+v' ) 


V*>  '  V°> 


ft’) 


-(l-2b-v‘ )/(b+v' ) 


+  A(0)61  j  I  ^  T/ 


q  ^  v'/(b+v')  .  v  -(l-2b-v')/(b+v') 


We  now  multiply  (67)  by  A^Aj  when  A^  is  an  arbitrary  vector  and  find 

_  /„  \  -l/(b+v') 


_  _ _  /  q 

A1u1uj(t)Aj  *  (A-ur(0)(  ^pTj  + 

2  r/«  \*2b/(b+v')  /q  \ 

•H(k')  -ft-) 


•l/(b+v' ) 


From  (57)  we  see  that 


t2-  t  >  1  (b  +  v’  >  0) 
o 
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Sufficient  for  the  left-hand  side  of  (69)  to  be  positive  is 

\  -2b/(b+v' )  f  v-l/(b+v') 

)  -(<■) 

which  requires,  using  (70), 

2b  <  1  (72) 


A  similar  analysis  applies  to  Aij  ;  however,  no  further  restrictions  on 
the  parameters  is  found. 

C.  An  Exact  Solution  of  the  Equations  with  Shear 
To  obtain  a  solution  of  the  equations  with  shear,  we  let 


V  e 


au't 


L  e 


au't 


^.w1e2lU't 
^•W2  e2aU,t 
*  W3  e2aU'1 


“l“2  *  W4  e 


2aU't 


A11  *  L1  e 
^22  *  L2  e 


au't 


au't 


A33  *  L3  e 


au't 


‘12 


L4  e 


au't 


(73) 


Substituting  these  forms  into  the  differential  equations  of  III. A,  we  find 
that  the  exponentials  cancel  and  that  an  algebraic  set  of  equations  for  the 
amplitudes  are  obtained.  It  is  possible,  with  some  algebra,  to  solve  the 
amplitude  equations  explicitly  In  terms  of  the  parameters  b  and  v'  . 

The  energy  components  are 


U1  _  1  +  6v'  +  4b 
7  Sir  ♦  2vr) 


(74) 


1  -  2b  1 

"TTTF 


(75) 


The  scale  components  are 


A11  _  1  +  6v'  +  4b 
“A - TTTP — 


(76) 


19 


*22  _  A33  _  1  -  2b 
~r  ~  l  +Tvr 

Me  see  that 

. 

A11  .  ^1_  ,  l±ivLlik 

h~2  — 5  1  -  2b 

The  off-diagonal  components  are 

~r  ‘ "  tt-2^-^3(1  ■ 2b){b  +  y,) 

The  Corrsin  parameter  is 


C  = 


ulu2 


~7  T 
U1  U2 


The  ratio  of  the  two  times  is 


1 

a 


q  .  1 _  /  1  -  2b 

ITA  "  1  +  2v'  / 3(b  +  v*) 


and  the  growth  rate,  a  ,  is 


a 


1  .  v*  /  1  -  2b 
v  o  1  V  2vr  /  3(b  +n 


We  notice  two  additional  Interesting  parameters: 


I  uiu2  I  u* 


a*Br  *  b  +  v' 


hi. 

A 


55. 


1  A 

J5 


i^ 


(77) 

(78) 

(79) 

(80) 

(81) 

(82) 

(83) 

(84) 

*5) 
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In  terms  of  the  anisotropy  parameter  o  defined  in  Ref.  2,  p.  4,  we  note 
that  this  condition  corresponds  exactly  to 

o*3 

The  value  of  o  that  has  given  satisfactory  results  for  a  wide  class  of 
turbulent  flows  is 


We  add  the  following  observations : 

1)  Several  numerical  integrations  suggest  that  any  solution  that  is  statisti¬ 
cally  realizable  initially  will  remain  so  and  will  asymptote  to  the 

exact  solutions  given. 

2)  If  the  scalar  scale  equation 

f?  *  c  T  V2  U‘  +  v'q  (c  0.35)  (86) 

is  adopted  instead  of  the  tensor  equation,  an  exponential  solution 
exists  and  has  the  same  qualitative  features  insofar  as  scale  and 
energy  are  concerned  (of  course,  no  scale  directivity  results).  This 
indicates  that  the  models  for  scale,  which  were  not  designed  to  fit 
homogeneous  shear  data,  are  quite  stable. 

With  only  two  adjustable  constants,  the  model  covers  qualitatively  both 
types  of  turbulent  flows.  We  have  shown  this  by  exhibiting  explicit  analy¬ 
tic  solutions  with  several  of  the  desired  features.  In  particular,  the 
analytic  solutions  for  the  homogeneous  turbulence  models  show  the  presence 
of  two  distinct  time  scales  which  characterize,  respectively,  the  rapid 
settling  of  the  tensor  character  of  the  flow  to  an  asymptotic  state  and  the 
slower  development  of  the  energy  and  mean  scale.  It  is  found  that  for  both 
grid  and  homogeneous  shear  turbulence,  the  ratio  of  the  two  scales  is  about 
ten. 

The  fast  time  is  the  redistribution  time,  A/q  ,  while  the  slower  one  is 
bA/q  (dissipation  scale)  for  grid  turbulence  and  v'A/q  (merging  scale) 
for  shear  flows  after  an  initial  transient.  This  feature  of  the  model 
solutions  seems  to  be  well  reflected  in  the  data. 
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V.  RELATION  BETWEEN  THE  HIGH  REYNOLDS  NUMBER  DISSIPATION 
AND  THE  EDDY  SIZE  REARRANGEMENT 


A  further  test  of  our  modeling  consists  In  comparing  the  spectra  predicted 
by  the  R-j j  model  with  experimental  data.  The  simplest  flow  on  which 
spectral  Information  (Energy  Spectra)  is  available  Is  grid  turbulence.  A 
first  cut  at  the  analysis  can  be  made  by  assuming  that  grid  turbulence  is 
Isotropic  (Ref.  5).  This  assumption  is  In  fact  not  correct,  and  a  more 
general  framework  is  needed  in  order  to  understand  the  full  details  of  the 
measurements.  A  remarkable  result  can  be  demonstrated,  however.  For  high 
turbulence  Reynolds  numbers  (qA/v  large),  the  energy  dissipation  is  given 
experimentally  by 


b  *  0.125 


(87) 


A  theorem  due  to  G.I.  Taylor  (Ref.  4,  p.  100)  states  on  the  other  hand 
that 

|£.-2v£  (88) 

The  two  statements  of  dissipation  are  equivalent  if  the  following  relation 
(Glushko,  Ref.  6)  is  assumed 


A 


2 


(89) 


We  demonstrate  below  that  a  simple  model  of  eddy  size  rearrangement  needed 
to  obtain  closure  of  the  R-j j  equation  implies  (89).  We  introduce  wave 
number  space  by 

*ij(*}  *  — y  JRij(?)  e1*’r  d*  (90) 


and  the  three-dimensional  spectrum,  E  ,  by  the  Karman-Howarth  relation 


5.  Cerasoli,  C.,  Donaldson,  C.duP.  and  Sandri,  G.,  "Fundamental  Research 
in  Turbulent  Modeling,"  AFOSR  TR-80-0324,  February  1980. 

6.  Glushko,  G.S.,  "Turbulent  Boundary  Layer  on  a  Flat  Plate  in  an 
Incompressible  Fluid,"  Bull.  Acad.  Sciences  USSR,  Mech.  Series  No.  4, 
1965,  pp.  13-23. 


The  normalization  is  chosen  so  that 


i-r 


E(k)  dk 


The  equation  becomes 

||  *  T(k)  -  2vk2E  (93) 

where  the  transfer  function  T  which  represents  "cascading"  (more  properly, 
eddy  size  rearrangement)  satisfies 


T(k)  dk  =  0 


i.e.,  eddy  size  rearrangement  does  not  change  the  energy.  A  simple  model 
for  T  ,  which  has  derivatives  only  (no  integrals)  both  in  k  space  and 
for  Rij  ,  is 


T  *  "  A  a[!  E  +  k  H]  +  I  A  °  ^ 


The  coefficients  have  been  chosen  so  that  for  Kolmogoroff  equilibrium, 
i.e.,  T  =  0  ,  the  only  solution  is 

E  =  const  •  k’5^3 


and,  in  addition,  Eq.  (94)  is  satisfied.  For  purposes  of  calculational 
feasibility,  it  is  most  desirable  that  the  model  equations  be  differential. 
If  both  and  R^j  are  to  be  governed  by  differential  equations,  then 

it  can  be  snown  that  all  terms  in  the  modeling  of  T  must  be  of  the  form 
apkn3nE/3kn  .  A  general  consequence  of  such  assumptions  is  that  $(k) 
will  exhibit  a  region  where  <j>  is  proportional  to  a  power  of  k  .  From 
the  point  of  view  of  theory,  it  seems  that  the  simplest  assumption  one  can 


the  point  of  view  of  theory,  it  seems  that  the  simplest  assumption  one  can 
make  on  the  eddy  size  rearrangement  model  Is  that  It  be  local  simultaneous¬ 
ly  in  wave  number  and  relative  separation.  To  obtain  Eq.  (89),  we 
consider  a  steady  source  at  k0 


Hi 


E  +  k 


+  2vk2E 


Q06*k 
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The  exact  solution  Is 


(q/A)ak 


e*  H 


ft--) 


-(k-kc)4K 


(k/k) 


57 T 


(97) 


with  H  the  Heaviside  function  and 


kft2A2  ^ 
*  *  “SRi - 3~~ 


(98) 


"We  now  impose,  as  integral  constraints,  the  definitions  of  q^/2  ,  X  and 
A  . 


.2  i 

r 

II 

|  E(k)  dk  =  E^IjU) 

(99) 

^  7 ‘J 

0 

f  A*"  E0k02 12(<) 

0 

(100) 

f  !<».  £„!,(«> 

(101) 

"  0 

where  we  have  Introduced  three  convenient  integrals. 

II  *  12  »  ^3  *  an<^ 

abbreviation 

r  _  ^o  _ < 

Eo  (q/A)  k".  * 

(102) 

After  some  manipulation,  we  find 


f?=rfjT^TaRe  (103) 

Numerical  evaluation  of  Ij  then  gives 

X2  *  —5 - ^ -  (104) 

(ir /4)  +  (oRe/3) 
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(105) 


Comparing  Eq.  (104)  with  (89),  we  obtain 

_2 

a  *  \  *  2.47  .  b-f 

It  is  interesting  to  note  that  the  empirically  accepted  values  of  a  lie 
between  3.25  and  2.5.  We  note  here  a  general  result,  namely  that  a  model 
of  the  eddy  size  rearrangement  term  in  the  R-j j  equation  determines  the 
asymptotic  model  in  the  UJuT  equation.  Thisdis  not  obvious  when  one 
merely  considers  the  contraction  of  the  R^j  equation  into  the  ujuj 
equation.  From  dimensional  arguments,  it  can  be  shown  that  a  family  of 
models  should  yield  the  Glushko  relation  and,  in  addition,  that  if  the 
spectrum  has  the  Kolmogoroff  form 

E - *  A  e2/3k_5/3  for  k  >  k 


then  the  dissipation  parameter  b  is  given  by 

b  *  -  "yy  A" 3/2  *  0.121A"3/Z  (106) 

5*3 

Experimentally,  A  •  1;  thus  Eq.  (106)  yields  a  value  for  b  that  is  close 
to  the  generally  accepted  value  of  b  ,  i.e.,  0.125. 

It  is  of  interest  to  observe  that  a  model  in  the  sense  of  the  familiar  Krook 
model  (or  single  relaxation  time  model)  would  not  require  any  derivatives 
at  all.  Unfortunately,  it  does  not  seem  possible  to  implement  such  a 
concept  in  modeling  cascading  because  the  "equilibrium"  spectrum  (assumed  of 
Kolmogoroff  form)  has  infinite  dissipation. 
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VI .  CONCLUSIONS 


We  have  tested  our  formulation  of  the  modeled  equations  for  the  R-j j  tensor 
by  Illustrating  Its  properties  In  the  case  of  homogeneous  turbulence.  Our 
analysis  of  the  model  of  homogeneous  turbulence  has  brought  out  four  main 
conclusions. 


1.  The  model  contains  only  two  adjustable  parameters  which  we  chose  to 
fix  from  grid  turbulence  data.  Good  qualitative  agreement  with  homo¬ 
geneous  shear  flow  results. 

2.  The  model  Implies  the  existence  of  two  distinct  time  scales  which  are 
separated  by  a  factor  of  about  10.  They  appear  after  an  Initial 
transient  phase  has  died  out  during  which  q/A_  becomes  approximately 
equal  to  U*  .  On  the  fast  scale,  A/q  z  (U ' )“ 1 ,  the  normalized 
deviator 


locks  Into  a  constant  value  Indicated  by  a  Corrsln  parameter  (Refs.  7 
and  8) 


-  0.5 


or  by  a  Bradshaw  number 


On  the  slower  scale,  v'(A/q)  -  0.07(A/q)  ,  the  energy  components  and 
the  mean  scale  grow  exponentially.  This  solution  can  be  thought  of  as 
a  superequilibrium  with  convection.  A  qualitative  physical  picture  Is 
as  follows.  During  the  Initial  transient,  when  (q0/A0)«U'  ,  the 
shear  brings  the  sudden  dlstortlon-llke  turbulence  up  to  convective 
equilibrium  (q/A  *  U'j,  while  when  (qo/Ao)  »  U'  ,  the  turbulence  decays 
grid-like  to  convective  equilibrium.  Then  a  merging  mechanism  takes 
over  (a  multi-layer  Brown-Roshko  effect)  so  that  eddies  fold  with  each 
other,  making  larger  ones  Indefinitely  (as  long  as  the  Imposed  shear 
provides  the  energy  to  sustain  the  process).  Once  the  merging  process 
takes  over,  the  eddy  structure  remains  fixed  and  exhibits  highly  direc¬ 
tional  Integral  scales  (see  4  below). 

7.  Harris,  V.G. ,  Graham,  J.A.H.  and  Corrsln,  S.,  "Further  Experiments  In 
Nearly  Homogeneous  Turbulent  Shear  Flow,"  J.  Fluid  Mech.  81,  1977,  pp. 
657-687.  Corrigendum,  J.  Fluid  Mech.  86,  1978,  pp.  795-795". 

8.  Corrsln,  S.  and  Kollman,  W.,  "Preliminary  Report  on  Suddenly  Sheared 
Cellular  Motion  as  a  Qualitative  Model  of  Homogeneous  Turbulent  Shear 
Flow,"  Proc.  SQUID  Symposium  on  Turbulence  In  Internal  Flows ,pp.  11-33 
(S.N.B.  Hurthy,  ed.),  Hemisphere  Publishing  Co.,  1977 
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3.  In  the  model,  the  two  transverse  energy  components  do  not  separate 
while  the  data  Indicate  that  such  a  separation  occurs.  The  splitting 
can,  however,  be  brought  Into  the  model  by  assuming  a  tensor  coeffi¬ 
cient  In  the  redistribution  equation.  We  have  taken  this  point  to 
constitute  a  refinement  at  this  stage  of  analysis. 


4.  The  calculated  angular  averaged  integral  scales  are  quite  directional. 
The  model  thus  gives  a  picture  of  the  eddy  structure  of  a  sheared 
turbulence.  With  our  choice  of  parameters 


0.19 


this  picture  could  eventually  be  tested  by  experiment. 


The  derivation  of  the  scale  equation,  the  agreement  of  the  model  with  homo 
genous  shear  flows,  and  the  derivation  of  the  Glushko  relation  are,  in  our 
opinion,  strong  arguments  in  support  of  developing  a  full  model  for  . 
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APPENDIX  A. 

THE  INTERCOMPONENT  REARRANGEMENT  MODEL  FOR  ISOTROPIC  TURBULENCE 


We  prove  that  our  modeling  of  the  two-point  intercomponent  rearrangement 
term  makes  this  term  vanish  automatically  for  isotropic  turbulence.  We  also 
show  why  this  statement  is  nontrivial  for  two-point  as  contrasted  to  one- 
point  models. 


The  intercomponent  rearrangement  two-point  tensor  is  given  by 

/  .  3u!(y)  3uj(x)  „ 

— p'(y) 


3xci 


3xcj 


A(J(xc,r)  - 

We  model  generally  Aij  as  (Eq.  (23)  in  the  text) 

A1jS-  j[RiJ-jVaatC1j] 


(A.l) 


(a. 2) 


From  the  definition,  Eq.  (A.l)  and  the  continuity  relations 


3u.(x)  Su^y) 


3xci  3xci 


=  0 


(A. 3) 


we  find  that  A^  must  be  traceless;  i.e., 

Akk  =  0 

From  Eq.  (A. 4)  we  conclude  that  C^  itself  must  be  traceless. 

*  J 


0  =  A, 


kk 


-*[ 


Rkk  *  Roc  +  Ckk 


;]  =  ‘  A  Ckk 


(A. 4) 
In  fact, 

(A. 5) 


For  homogeneous  turbulence,  R-jj  ,  C-|  i  and  hence  A^j  depend  on  the^. 
separation  r  only.  We  Fourier  transform  Eq.  (A. £)  with  respect  to  r 
using  Eq.  (90)  and  analogous  equations  to  define  Aii  and  Cii  .  We  find 
in  wave  number  space 

VR>  ’  -?[♦»/>- 5  Vaa  +  CtJ(k)]  (*•«) 


For  all  the  models  that  we  consider  for  Isotropic  turbulence,  S-u  satis¬ 
fies  continuity  automatically.  We  can  then  determine  C^j  uniquely.  In 
fact,  from  Isotropy  the  most  general  symmetric  second  rank  tensor  can  be 
written  In  terms  of  two  scalars,  a(k)  and  b(k)  as 

C^dT)  -  aflcjdjj  ♦  btk^kj  (A. 7) 
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where  £  denotes  the  unit  vector 

kj  *  k./k  (A. 8) 

But  C^j  Is  traceless  by  Eq.  (A. 5);  thus,  from  Eq.  (A. 7) 

0  -  3a  +  b  (A. 9) 

Substituting  (A. 9)  Into  (A. 2),  we  have 


(A. 10) 

Continuity  of  A^j  requires 

0  .  k,AfJ 

■  '  A  k1  [*1j(k)  ■  7  stj»ao  +  Sj] 

‘  -  *  kj  [-  5  ♦.  •  2k<k>] 

(A. 11) 

Hence  a(k)  is  determined: 

•<k>  is 

(A. 12) 

and  from  Eq.  (A. 10),  Is  determined  uniquely  by  the  energy  spectrum 

We  now  substitute  (A. 13)  Into  (A. 6)  and  use  the  Karman-Howarth  theorem, 

Eq.  (91)  In  the  form 

♦lj(k)  *  (61j  '  Vj)  "P  (A-14) 

We  find 

A1j(k)  *  "  ?[(61j‘C1kj  )  “T"!  61j*oa  -(61j‘3k1kj)  V] 

*  0  (A. 15) 
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This  completes  the  proof  of  our  theorem. 


If  we  had  assumed  a  model  for  Ajj  without  the  term,  i.e., 

A1j  '  -  A  [  RU  -  T  <*•“> 

we  would  have  violated  tlje  requirement  that  A^j  vanish  for  isotropic 
turbulence  because  (in  r  space),  using  the  Karman-Howarth  representation 


R1j  ■  4  [?Sij  +  <f  ■ 


(A. 17) 


We  find 


Ru 


-  4 

3  ij  aa 


=  18  (6ij  "  3?i?j)r  f’(r)  *  0  (A*18) 


Following  von  Karman  and  Howarth,  we  have  considered  here  only  mirror 
symmetric  isotropic  turbulence  and  hence  neglected  the  skew-symmetric  part 
of  R^j  (and  of  C-jj  ).  A  straightforward  calculation,  which  follows 
essentially  the  steps  given  here,  shows  that  the  skew-symmetric  parts  of 
♦ij  an<l  *ij 


pij 

\/ 


Ua 


£  A(k) 


“ij 

s/ 


Mjo1 


k  c(k) 


(A. 19) 


are  related  by 


A(k)  +  c(k)  =  0  (A. 20) 

in  order  to  satisfy  the  theorem.  The  result  is  therefore  valid  for  general 
isotropic  turbulence. 

Mathematically  we  may  summarize  our  analysis  in  the  following  theorem: 

Any  isotropic  (second  rank)  tensor  that  is  traceless  and  divergence- 
free  Is  skew- symmetric  (and  therefore  a  curl). 

Even  though  it  Is  possible  that  alternative  models  satisfy  the  Isotropy 
requirement,  i.e.,  A^<  *  0  ,  It  Is  an  advantage  of  the  model  proposed  here 
that  it  does  so  explicitly. 
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I.  INTRODUCTION 


The  study  of  turbulent  flows  began  over  a  century  ago  with  the  work  of 
Reynolds.  Since  that  time,  our  understanding  of  turbulent  flows  has 
grown,  but  not  as  rapidly  or  dramatically  as  In  many  other  fields  of 
science.  A  very  old  and  useful  concept  In  the  study  of  turbulence  Is 
the  turbulent  length  scale,  or,  more  physically,  the  "eddy"  size.  The 
notion  of  a  typical  eddy  size  Is  Intimately  connected  with  the  turbulent 
transport  of  all  flow  variables,  and  Is  an  aid  In  qualitatively  under¬ 
standing  turbulent  flows.  The  concept  of  length  scale  enters  explicitly 
when  one  attempts  to  simulate  turbulent  flows  numerically  by  use  of 
"second  order  closure"  models.  Data  on  the  magnitude  of  the  turbulent 
scales  Is  required  as  Input  for  such  modeling  schemes,  and  the  present 
experimental  study  provides  data  on  turbulent  scales  for  an  annular 
flow  field. 

The  manner  In  which  length  scales  appear  in  closure  models  will  be 
discussed  in  Section  II.  We  will  also  present  the  definition  of  the 
correlation  tensor  and  the  way  in  which  length  scales  can  be  defined 
via  this  tensor.  The  experimental  facilities  and  techniques  used  here 
will  be  discussed  in  Section  III,  while  data  on  the  basic  flow  field  (mean 
and  turbulent  velocities)  will  be  presented  in  Section  IV.  Section  IV  will 
also  include  a  discussion  of  theoretical  concepts  related  to  the  mean  and 
turbulent  velocity  fields,  and  a  comparison  of  our  data  with  that  from 
previous  experiments.  Turbulent  length  scales  obtained  by  use  of  auto¬ 
correlation  measurements  and  Taylor's  hypothesis  are  presented  in  Section  V, 
where,  again,  theory  and  experiment  are  discussed.  The  notion  of  the 
von  Karmon  length  scale  is  also  introduced  and  compared  to  data. 

Section  VI  contains  results  obtained  from  two-point  correlation  measure¬ 
ments,  and  data  are  compared  to  previous  experiments  and  our  theoretical 
knowledge.  In  particular,  results  are  discussed  using  concepts  developed 
for  Isotropic  and  homogeneous  flows.  Although  the  annular  flow  field  is 
neither  isotropic  nor  homogeneous  globally,  these  concepts  are  useful 
and  consistent  with  our  data  In  certain  spatial  regions  where  the  mean 
flow  gradients  are  small. 

The  two-point  correlation  equations  are  also  presented,  and  we  qualitatively 
discuss  the  effect  of  the  production  terms  on  the  various  correlation 
functions.  The  results  from  this  study  are  summarized  in  Section  VII 
and  future  research  topics  are  discussed. 
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II.  SECOND  ORDER  CLOSURE  MODELS 
AND  THE  CORRELATION  TENSOR 


The  typical  method  for  solving  turbulent  flows  requires  the  decomposition 
of  dependent  flow  variables  into  mean  and  fluctuating  parts.  The  total 
velocity  field  is  written  as 


where  is  the  mean  velocity  (either  a  time  or  ensemble  average)  and 

u.j  is  the  fluctuating  part;  the  subscript  i  runs  from  1  to  3  for  the 
three  velocity  components.  A  similar  separation  is  used  for  all  other 
dependent  variables,  such  as  temperature,  density,  etc.  The  following 
equations  result  when  such  a  decomposition  is  applied  to  the  Navier- 
Stokes  equations  for  incompressible,  constant  temperature  flow: 


3Uj  3U. 
-3F+UJ  aT  + 

J 


1 

P  3xi 


-  W2Ui 


where  P  and  p  are  the  mean  pressure  and  density,  respectively,  v 
is  the  kinematic  viscosity,  repeated  indices  are  summed  over  and  the 
overbar  denotes  an  average.  We  find  that  the  equation  for  mean  variables 
involves  the  correlation  of  fluctuating  variables,  u-jUj  .  The  equation 
for  u^uj  reads 


Here  we  find  that  the  equation  for  the  second  order  quantity,  u^u.  ,  involves 
third  order  variables,  UjUju^  ,  pressure-vel oci ty  correl ati ons ,  J 
puj  ,  and  pressure-velocityJgrad1ent  terms,  p(  u^/  xj)  .  One  can  continue 
this  process  and  the  equation  for  third  order  variables  will  Involve  fourth 
order  terms.  This  problem  of  always  having  one  equation  less  than  the 
number  of  dependent  variables  Is  the  "closure"  problem.  The  equations 
can  be  closed  at  the  first  order  by  writing  u j uj  as  a  function  of  the 
mean,  or  first  order,  variables.  This  is  usually  done  via  the  concept  of 
an  eddy  viscosity,  vE  ,  such  that 


This  method  essentially  treats  turbulent  mixing  in  a  manner  similar  to 
viscous  diffusion  and  can  yield  reasonable  results  for  certain  flow  fields 
The  next  level  of  sophistication  is  to  close  the  system  at  second  order, 
and  this  allows  the  physics  of  turbulence  to  be  more  accurately  modeled 
than  the  use  of  eddy  viscosities.  The  set  of  equations  are  closed  at 
second  order  by  writing  equations  for  u-f ujuk  ,  pui  and  p(3uj/3xj) 
solely  in  terms  of  second  order  quantities,  ihe  present  A.R.A.P.  model 
uses  the  following  equations: 
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where  the  constants  Vc  and  Pc  are  determined  experimentally.  The 
symbol  A  represents  a  variable  length  which  is  locally  representative 
of  the  Integral  turbulent  length  scale,  or  alternately,  a  typical  eddy 
size. 


The  determination  of  A  requires  knowledge  of  the  correlation  tensor, 
defined  as 


Rij  (r.;  Ar) 


=  ui (r )  Uj (r  +  Ar ) 


For  the  case  Ar  s  0  ,  R-j j  (r;  0)  equals  the  turbulent  intensity  squared 
for  each  velocity  component  when  i  =  j  ,  and  the  turbulent  shear  stress 
when  i  i  j  .  R-j  j  (j^,  AjO  will  in  general  be  a  function  of  the  vector 
position  r  within  the  flow  field  and  the  vector  separation  Ar  .  A 
normal izecTvers ion  of  ftjj  will  be  used  in  most  of  this  work,  defined  as 


Ri j (n;  Ar)  =  Rij(r;  Ar) 


Correlation  lengths  may  be  generated  from  R-j  j  in  a  number  of  ways. 
One  method  is  to  first  perform  an  angular  average  over  all  separation 
directions;  this  is  known  as  a  "Batchelor  average"  and  is  defined  by 


R-  Ar}  Ar)  dft 

*  J 


where  dft  is  the  solid  angle  vector.  A  length  scale  may  be  defined  by 
integrating  over  the  scalar  separation  variable  Ar  , 


L 


r 


Ar 


d(Ar) 


Note  that  one  may  generate  nine  independent  length  scales  which  are 
functions  of  position  r  ,  and  that  Li j  is  a  tensor.  Another  way  to 
define  length  scales  is  to  integrate  over  a  single  separation  direction 
while  holding  the  other  two  separations  at  zero.  We  take  a  case  using 
cylindrical  coordinates  as  an  example  where 


and 


r  =  (r,  <p,  z) 


r  =  (Ar,  A<t>,  Az) 


Let  Ar  *  A<J>  =  0  ,  and  a  length  scale  A^,  z  can  be  defined  as 
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R^r;  o ,o ,Az )  d(Az) 


where  the  subscripts  indicate  the  tensor  component  and  the  integration 
direction,  respectively.  Note  that  is  not  a  tensor  and  the 

subscripts  are  used  purely  for  bookkeeping  purposes.  A  total  of  27 
length  scales  can  be  defined  by  the  preceeding  technique  (9  components 
times  3  directions).  Second  order  closure  models  typically  require  a 
single  length  scale,  and  the  richness  of  information  contained  in  R^j 
and  its  associated  scales  must  be  replaced  by  a  single  number.  This 
is  best  done  when  one  understands  the  structure  and  details  of  the 
correlation  tensor. 

The  present  study  involves  axisymmetric,  annular  flow  and  cylindrical 
coordinates  are  used,  where  (r,<f>,z)  are  the  radial,  azimuthal  and 
axial  coordinates.  The  corresponding  velocity  components  are 
(u,v,W  +  w),  where  we  have  used  the  fact  that  the  mean  radial  and 
azimuthal  velocities  are  negligible  and  taken  as  zero  (U  =  V  =  0). 

The  flow  is  assumed  to  be  well  developed  and  axisymmetric,  so  R-jj 
is  independent  of  <p  and  z  ,  and  we  write  R. .  as 

*  J 


Rjj  (r;  Ar,A<j> ,Az) 

Azimuthal  separations  were  held  to  zero  in  all  experiments,  and  either 
Ar  =  0  ,  Az  i  0  ,  or  Ar  f  0,  Az  =  0  ,  so  that  the  separations  were 
either  pure  axial  or  pure  radial.  The  axial  separations  were  implied 
by  use  of  Taylor's  hypothesis  which  will  be  subsequently  discussed. 

The  axisymmetric  nature  of  this  flow  results  in 


V-jUg  =  u-|V2  =  v-|W2  =  v2w-|  =  0 

where  the  subscripts  refer  to  position, 

1  =  r 

2  =  r  +  Ar 

We  therefore  have  a  tensor  with  5  independent,  non-zero  components  for  the 
annular  flow  field. 
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III.  EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 


The  wind  tunnel  used  in  this  study  is  described  in  a  previous  A.R.A.P. 

Report  (Ref.  1),  and  a  brief  description  of  the  facility  will  be  given 
here  for  completeness.  The  facility  is  shown  schematically  in  Figure  (1), 
where  the  test  section  is  an  annulus  with  Inner  and  outer  radii ,  R»  and 
Ro,  of  2.54  and  21.59  cm,  respectively.  The  distance  from  the  end*of  the 
contraction  section  to  the  test  section  can  be  varied  by  the  addition  or  sub¬ 
traction  of  cylindrical  sections,  and  the  minimum  and  maximum  distances  can 
be  varied  from  15  to  45  diameters  (2Ro) -  The  contraction  section  admits  air 
through  the  region  shown  in  Figure  (1)  and  vanes  are  situated  to  direct  the 
incoming  air.  The  vanes  may  be  set  at  any  angle  from  0°  to  45°,  and  when 
the  vane  angle  Is  non-zero,  angular  momentum  Is  Imparted  to  the  incoming  air. 

The  contraction  process  will  then  create  a  swirling  flow  in  the  annulus.  A 
constant  speed  motor  was  used  In  conjunction  with  a  variable  speed  fan  drive 
to  provide  flow  through  the  tunnel,  and  the  mean  axial  velocity  could  be 
varied  between  approximately  1200  and  3000  cm/sec.  Before  discussing  the 
data  acquisition  techniques,  we  will  discuss  certain  problems  associated  with 
the  wind  tunnel. 

Initial  measurements  made  by  this  author  showed  that  the  flow  field  was  highly 
unsteady.  This  was  apparent  In  the  behavior  of  the  turbulent  velocity  power 
spectra  as  low  frequencies  were  approached  (f  •*  0).  Typically,  as  f  -*•  0,  the 
power  spectral  density  approaches  a  constant  value  which  can  be  used  to  define 
a  turbulent  length  scale.  Our  Initial  experiments  showed  a  divergent  behavior 
where  the  spectral  level  Increased  as  f  -*  0,  and  this  Indicated  a  high  degree 
of  unsteadiness. 

The  source  of  unsteadiness  was  traced  to  the  Inlet  section.  Initially,  a  fine 
screen  was  placed  around  the  inlet  *nd  did  not  provide  any  damping  of  atmos¬ 
pheric  motions.  Such  motions  woulo  oe  amplified  as  the  flow  contracted  and 
gave  rise  to  the  unsteadiness;  the  problem  was  remedied  by  placing  a  5  cm 
thick  piece  of  polyurethane  foam  around  the  Inlet  section.  This  provided  a 
pressure  drop  across  the  foam  of  more  than  10  times  the  dynamic  pressure  based 
on  the  inlet  velocity  (pU|/2,  p  Is  the  fluid  density  and  Uj  Is  the  inlet 
velocity),  and  filtered  out  all  the  atmospheric  motions. 

A  detailed  survey  of  the  mean  and  turbulent  velocity  fields  was  taken  following 
the  solution  of  the  unsteadiness  problem.  These  measurements  revealed  that 
the  flow  field  was  not  axl symetrl c .  One  Indication  of  this  lack  of  symmetry 
was  that  the  maximum  mean  axial  velocity  was  at  r  ■  .55Ro  on  the  near  side  of 
the  Inner  tube  U  *  0°),  while  the  maximum  mean  occurred  at  r  »  .35R«  on  the 
far  side  of  the  Inner  tube  ($  *  180° ) .  This  behavior  was  believed  to  be  asso¬ 
ciated  with  the  shape  of  the  contraction  section,  which  may  have  caused  the  flow 


1.  Bllanln,  Alan  J.,  Snedeker,  Richard  S.,  Sullivan,  Roger  D. ,  and  Donaldson, 
Coleman  duP.,  "Final  Report  on  an  Experimental  and  Theoretical  Study  of 
Aircraft  Vortices,"  AF0SR-TR- 75-0664,  A.R.A.P.  Report  No.  238,  February 
1975. 


to  become  separated.  The  problem  was  remedied  by  placing  three  sections  of 
honeycomb  after  the  contraction  section,  as  shown  in  Figure  (2).  The  pressure 
drop  associated  with  the  honeycomb  was  great  enough  to  reduce  the  asymmetries 
created  in  the  contraction  section,  and  data  showed  that  the  flow  was  symmetric 
to  within  our  resolving  capabilities.  The  disadvantage  of  this  solution  was 
that  swirl  cannot  be  introduced  by  the  vanes.  Swirl  now  must  be  introduced 
after  the  honeycomb  In  the  annulus  region,  or  the  contraction  section  must  be 
redesigned  to  eliminate  the  asymmetries.  The  solution  of  using  honeycomb  was 
the  quickest  and  most  Inexpensive  one,  and  enabled  us  to  make  a  very  complete 
set  of  measurements  on  non-swirling  annular  flow. 

Fluid  velocities  were  measured  using  hot  film  X  probes  which  allowed 
parallel  and  traverse  velocities  to  be  determined.  The  probes  were  manufac¬ 
tured  by  Thermo  Systems  Incorporated  (TSI,  model  number  1240-20),  and  the 
sensor  elements  were  .05  mm  In  diameter  and  1  mm  long.  TSI  anemometers 
and  llnearlzers  (model  1054)  were  used  In  conjunction  with  the  X  probes. 

A  variety  of  probe  holders  were  designed  and  used  In  this  study,  but  only  two 
need  to  be  described.  The  first  probe  holder  allowed  a  single  probe  to  be 
traversed  radially  on  either  side  of  the  Inner  tube.  This  was  used  to  obtain 
mean  and  turbulent  velocities  as  a  function  of  radius  and  determine  the  sym¬ 
metry  of  the  flow  field.  A  TSI  sum  and  difference  amplifier  (model  1063) 
was  used  during  the  single  probe  measurements.  This  provided  output  voltages 
which  were  In  proportion  to  the  parallel  (sum)  and  transverse  (difference) 
velocities.  The  second  probe  holder  allowed  one  probe  to  remain  at  the  fixed 
radius,  rp,  while  a  second  probe  was  traversed  to  other  radii.  This  allowed 
two-point  correlations  to  be  measured  for  the  case  of  pure  radial  separations. 
No  sum  and  difference  amplifiers  were  used  In  the  two  probe  measurements  as 
the  signals  were  added  and  subtracted  In  the  data  analysis  programs.  The 
signals  from  the  anemometer  and  llnearlzer  units  were  amplified  by  AC  coupled 
amplifiers  with  a  bandwidth  from  2Hz  to  over  20,000Hz.  These  amplified  signals 
were  transmitted  via  cable  to  the  computer  facility  and  Into  a  low  pass 
filter. 

The  filtered  output  was  digitized  by  use  of  an  analog  to  digital  (A/D) 
converter.  Digitized  data  was  then  stored  on  tape  and  subsequently  analyzed. 
The  sample  rate  for  the  single  probe  experiments  was  5,000  samples  a  second, 
which  gave  a  Nyqulst  frequency  of  2,500Hz.  The  cutoff  frequency  for  the 
low  pass  filter  was  2,000Hz  and  this  eliminated  any  sorce  of  aliasing  error. 

The  two  probe  measurements  used  the  A/D  converter  In  the  multiplex  mode.  The 
sample  rate  for  all  four  data  channels  was  50,000  samples  per  second,  which 
yielded  a  sample  rate  per  channel  of  12,500Hz  and  a  Nyqulst  frequency  of  6250Hz 
the  low  pass  cutoff  frequency  was  5,000Hz.  The  data  aqulsltlon  system  Is 
shown  schematically  In  Figure  (3),  where  the  flow  of  data  from  probes  and 
anemometers  through  amplifiers  and  filters,  into  the  A/D  converter  and  eventu¬ 
ally  to  storage  on  tape  can  be  seen.  Computer  programs  were  used  to  calculate 
mean  velocities,  RMS  fluctuations,  autocorrelations,  power  spectra  and  corre¬ 
lation  functions. 

All  experiments  were  conducted  with  a  maximum  axial  velocity  of  1960  cm/sec 
*  and  a  mass  flow  velocity,  Wr,  (defined  In  Section  IV)  of  1660  cm/sec.  The 
flow  Reynolds  based  on  mass  flow  velocity  and  the  hydraulic  diameter, 
dh  •  2(R  -  Rj),  was 
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R  -  -  4.2  x  105 

e  v 

while  the  distance  from  the  last  section  of  honeycomb  to  the  test  section,  z, 

was 

z  »  30  D0  •  34  dh 


Three  sets  of  measurements  were  made  for  the  two  point  correlations,  where  the 
fixed  probe  positions  were  rc  *  .27,  .43  and  .75  .  (Note  that  from  this  point, 
all  distances  will  be  normalized  by  Rq  .)  These  three  positions  were  chosen 
because  the  mean  axial  velocity  gradient  (dW/dr)  was  positive,  zero  and  negative, 
respectively,  at  these  positions.  Such  measurements  enabled  us  to  examine  the 
effect  of  mean  shear  and  geometry  on  the  correlation  functions. 


IV .  MEAN  AND  TURBULENT  VELOCITY  MEASUREMENTS 


A  series  of  single  probe  measurements  were  made  to  map  out  the  mean  and 
fluctuating  velocities  In  the  annulus.  The  following  variables  were  measured 
as  a  function  of  radial  position:  the  mean  axial  or  downstream  velocity, 

W,  the  turbulent  Intensities  for  the  radial,  azimuthal,  and  axial  velocities, 
u,  v,  w  [*  (u7)'/2t  (7?)l/2,  (» J7)1/2]t  and  the  turbulent  shear  stress,  uw. 


Mean  Velocity 


Figure  4  shows  the  mean  axial  velocity  versus  radius  where  velocities  have 
been  normalized  by  the  mass  flow,  or  bulk,  velocity,  Wg,  defined  as 


W 


B 


W(r)rdr  dtp 


Rj  and  R0  are  the  Inner  and  outer  radii,  respectively.  Measurements  were 
made  on  each  side  of  the  center  tube  (denoted  as  "near"  and  "far"  side)  and 
W  Is  symmetric  to  within  the  accuracy  of  the  measurements.  It  Is  well  known 
that  the  velocity  can  be  represented  logarithmically  for  the  outer  region  In 
pipe  flow  (distance  away  from  wall  greater  than  .05Ro),  and  a  similar  behavior 
exists  for  the  present  case.  Figure  S  shows  the  axial  velocity  versus  dis¬ 
tance  from  wall  In  semi-log  coordinates,  where  the  upper  and  lower  graphs 
are  for  distances  away  from  the  Inner  and  outer  radial  walls,  respectively. 
Both  cases  show  good  logarithmic  fits  until  the  central  annular  region  Is 
approached;  the  velocities  fall  under  the  log  profile  for  distances  from  the 
Inner  radial  wall  and  are  above  the  log  profile  for  distances  from  the  outer 
radial  wall.  No  analytical  solutions  exist  for  this  annular  flow,  therefore 
no  comparisons  can  be  made  between  theory  and  experiment.  In  the  pipe  flow 
case,  the  log  profile  Is  an  excellent  fit  right  to  r  *  0,  the  zero  mean 
gradient  position  (Nlkuradse,  Ref.  2).  The  present  case  shows  deviations 
from  the  log  profile  near  the  zero  gradient  position,  r  *  .43,  but  one  must 
be  cautious  when  relating  the  behavior  at  r  ■  0,  dW/dr  *  0  In  pipe  flow  to 

that  at  r  t  0,  dW/dr  *  0  In  annular  flow.  The  point  r  ■  0  Is  very 

different  geometrically  from  a  point  such  as  r  *  .43. 


B.  Turbulent  Intensities 


The  root  mean  square  turbulent  velocities  are  presented  In  Figure  6  as  a 
function  of  radius.  The  velocities  were  normalized  by  the  bulk  velocity  and 
the  data  are  from  near  side  measurements.  As  In  the  case  for  the  mean  axial 
velocity,  near  and  far  side  turbulent  Intensities  agreed  to  within  measurement 
accuracy.  A  number  of  features  can  be  discussed  In  terms  of  our  general 


2.  Nlkuradse,  J.,  "Gesetzmablgkelt  der  Turbulenten  Stromung  In  Glatten 
Rohren,"  Forsch.  Arb.  Ing.-Wes.  No.  356,  1932. 


knowledge  of  turbulent  shear  flows.  In  the  regions  away  from  the  zero 
gradient  position  (r  *  .43),  we  have  w  >  v  =  u.  This  Is  a  common  feature 
where  the  velocity  component  receiving  energy  by  production  directly  (w) 

Is  greater  than  the  other  two  components.  The  velocity  component  perpen¬ 
dicular  to  the  gradient  direction  (v)  Is  greater  than  the  component  parallel 
to  the  gradient  direction  (u);  similar  behavior  Is  observed  and  discussed  by 
Champagne  et.al.  (Ref.  3).  The  parabolic  behavior  of  w(r)  Is  consistent 
with  pipe  flow  measurements  and  w  Is  a  minimum  at  the  zero  gradient  posi¬ 
tion  where  no  local  production  occurs.  Two  features  of  the  data  differ  from 
expectations.  The  first  Is  the  "bulge"  In  the  w  profile  near  r  -  .65. 

This  Is  a  reproducible  feature  and  we  can  offer  no  explanation  for  It  as  no 
anamolous  behavior  occurs  In  other  flow  variables  near  this  position.  The 
second  feature  Is  the  lack  of  Isotropy  at  the  zero  gradient  position,  that 
Is,  w  f  u  *  v  at  r  »  .43,  and  we  define  an  Isotropy  factor  y  =  (u/w)r«  43 
*  .72.  Note  that  the  measured  values  for  the  transverse  velocities  (u  and’v) 
were  corrected  for  the  effects  of  tangential  cooling  (see  Appendix  A).  Had 
this  not  been  done,  the  Isotropy  factor  would  have  been  even  lower.  The  degree 
of  tangential  cooling  Is  embodied  In  the  k  factor,  which  was  measured  for 
the  hot  film  probes  used  here.  The  values  agree  with  the  published  results  of 
Frleke  and  Schwarz  (Ref.  4)  and  Jorgensen  (Ref.  5),  but  a  fair  amount  of 
scatter  exists  In  all  the  measurements.  If  we  take  the  largest  k  measured 
(which  results  In  the  greatest  correction  factor),  the  Isotropy  factor  Is 
Increased  to  about  .78,  as  opposed  to  .72.  This  lack  of  Isotropy  may  be  a 
consequence  of  the  flow  field  not  being  fully  developed.  Priest  (Ref.  6) 
shows  data  on  the  evolution  of  w  and  u  downstream  In  pipe  flow,  and  there 
Is  a  tendency  for  w  to  become  Independent  of  downstream  distance  before 
u  (and  v).  Priest's  results  show  that  at  r  ■  0, 


.70 

z/D  »  20 

.76 

z/D  -  30 

.78 

z/D  >  45 

The  present  result  concerning  Isotropy  can  be  compared  to  pipe  flow  measure 
ments  at  r  ■  0,  and  the  values  of  y  are  summarized  below. 


3.  Champagne,  F.  H.,  Harris,  V.  G.  and  Corrsln,  S. ,  "Experiments  on  Nearly 
Homogeneous  Turbulent  Shear  Flow,"  J.  Fluid  Mech.  41_,  pp.  81-139. 

4.  Frleke,  C.  A.  and  Schwarz,  W.  H.,  "Deviations  from  the  Cosine  Law  for 
Yawed  Cylindrical  Anemometer  Sensors,”  J.  Applied  Mech.,  Dec.  1968, 
pp.  655-662. 

5.  Jorgensen,  F.  E.,  "Directional  Sensitivity  of  Wire  and  Fiber-film  Probes," 
Dlsa  Report  No.  U,  May,  1971,  pp.  31-37. 

6.  Priest,  A.  J.,  "Incompressible  Turbulent  Flow  In  Pipes  and  Conical  Diffusers, 
Ph.D  Thesis,  University  of  Salford,  England,  1975. 
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Annular  Flow  r«  .43  Pipe  Flow«0 

f . . .  . . \  t  —  -  ' 

Present  Experiment  Laufer(Ref.  7)  Priest(Ref.  6)  Sabot  and  Comte-Bel lot (Ref .  8) 


.72 


.98 


.78 


.90 


Typical  error  bounds  on  y  are  about  .04,  so  that  substantial  disagreement 
exists  In  the  literature  concerning  the  degree  of  anisotropy  at  pipe  center. 

All  three  pipe  flow  measurements  were  made  at  downstream  distances  greater  than 
90  pipe  diameters,  and  lack  of  fully  developed  flow  Is  not  responsible  for  the 
differences.  It  Is  not  clear  whether  the  various  authors  corrected  the  data 
for  tangential  cooling,  and  If  this  Is  the  cause  for  the  differences.  The  reader 
should  note  that  a  similar  problem  exists  In  circular  jet  measurements. 

Vfygnanskl  and  Fiedler  (Ref.  9)  found  y  ■  *86  on  jet  centerline  (and  does 
compensate  for  tangential  cooling),  while  Gibson  (Ref.  10)  found  perfect  Iso¬ 
tropy,  y  ■  1.  on  centerline. 

C.  Turbulent  Shear  Stress 

Measurements  of  the  turbulent  shear  stress,  uw,  and  the  shear  stress  coefficient, 
uw/(u?)'/2(w?)'/s  were  made  as  a  function  of  radius.  These  measurements  are 
sensitive  to  both  probe  calibration  and  orientation,  and  this  can  be  seen  by 
noting  that 


"  '  (eA  *  *B,(eA  4  «B>  ‘  -  'I 

where  e«  and  eg  are  the  fluctuating  voltages  from  sensors  A  and  B  of 
an  X  wire  probe.  If  (1)  the  sensitivities  of  the  two  sensors  are  not  equal, 
or  (2)  the  probe  axis  Is  not  parallel  to  the  mean  flow  direction,  the  measured 
value  of  uw  will  be  In  error.  Such  errors  can  be  substantially  reduced  by 
appropriately  averaging  two  measurements  where  the  probe  Is  rotated  180°  about 
Its  axis,  thereby  Interchanging  the  positions  of  sensors  A  and  B.  This 
procedure  was  used  and  results  for  the  shear  stress  coefficient  are  presented 
In  Figure  7.  Measurements  from  the  near  and  far  sides  are  shown  and  the  degree 


7.  Laufer,  J.,  "The  structure  of  Turbulence  In  Fully  Developed  Pipe  Flow," 
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8.  Sabot,  J.  and  Comte-Bellot,  G.,  "Intemlttency  of  Coherent  Structures  in 
the  Core  Region  of  Fully  Developed  Turbulent  Pipe  Flow,"  J.  Fluid  Mech.  74, 
1976,  pp.  767-96. 

9.  Hygnanskl,  I.  and  Fiedler,  H.,  "Some  Measurements  In  the  Self-Preserving 
Jet,"  J.  Fluid  Mech.  38,  1969,  pp.  577-612. 

10.  Gibson,  M.  M.,  "Spectra  of  Turbulence  In  a  Round  Jet,"  J.  Fluid  Mech.  15_, 
1963,  pp.  161-73. 
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of  asymmetry  Is  the  sane  order  as  the  experimental  uncertainties.  The  values 
of  approximately  0.4  In  the  constant  coefficient  regions  are  In  good  agreement 
with  the  value  of  0.4  given  by  Laufer  for  pipe  flow.  The  shear  stress  coeffi¬ 
cient  Is  zero  at  r  -  .43,  where  the  mean  gradient  Is  zero,  and  the  sign  change 
In  the  coefficient  relates  to  the  sign  change  In  dW/dr.  The  shear  stress  Is 
presented  In  Figure  8,  where  we  have  averaged  the  near  and  farjlde  data.  The 
dashed  line  represents  a  fit  to  a  theoretical  expression  for  uw  which  we 
will  now  discuss. 

The  derivation  for  uw  as  a  function  of  radius  In  pipe  flow  Is  a  textbook 
exercise  which  can  be  found  In  Fluid  Dynamics  by  Dally  and  Harleman  (Ref.  11). 
The  derivation  requires  the  assumption  of  well  developed  flow,  and  when 
applied  to  the  case  of  annular  flow  gives 
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u«,i  and  u*  o  *re  friction  velocities  based  on  the  mean  gradient  at  the  Inner 
and  outer  radii.  The  variable  r  In  the  above  expressions  Is  the  dimensional 
radius,  as  this  allows  one  to  see  how  Rq  and  Rj  appear  in  the  equations. 

The  case  Rj  *  0  yields 

A  ■ 


B  ■ 


and 


uw  ■  u2 
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dW 

Hr 


which  Is  the  well  known  result  for  pipe  flow.  The  friction  velocities  can  be 
related  to  the  downstream  pressure  gradient  by 


11.  Dally,  J.  H.  and  Harleman,  D.  R.  F.,  Fluid  Dynamics,  Addlson-Mesley 
Publishing  Co.,  1966. 
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and  the  expression  reduces  to  the  pipe  flow  result  for  Rj  -  0. 


The  theoretical  expression  for  uw  Involves  two  constants,  u*  0  and  u*  j, 
which  can  be  obtained  by  direct  measurement.  Boundary  layer  m&asuremenis 
were  not  made  to  determine  gradients  at  the  walls  and,  as  a  consequence,  the 
following  methodology  will  be  considered.  The  position  at  which  uw  ■  0 
can  be  used  to  determine  the  ratio  of  u*  0  to  u*  r  ,  and  using  uw  *  0 
at  r  *  .43,  we  find 


1.33 


uw  -  4.0  ['  -227  r  -  *227  &J  ♦  v  $ 

The  experimental  data  can  be  fit  to  the  above  expression  and  yields  u*>0  »  .033  UR. 
Alternately,  the  friction  velocity  can  be  computed  from  measured  values  of 
dP/dz  (again  using  u*  r  «  1.33  u*  0)»  which  gives  u*  0  *  .040  Wg.  Finally, 
the  data  of  Koch  and  Feind  (Ref.  12  J  show  u*>0  «  .039  An  for  annuli  with 
R]/Ro  *  -6  and  .8.  These  three  values  are  summarized  below. 


u*.o/“B 


Technique 


.033 

.040 


Fit  Shear  Stress  Data 
to  Theory 

Pressure  Gradient  Measurement 


Present 

Study 


.039  Koch  and  Feind 


The  fit  to  data  yields  a  u*  0  below  that  of  the  other  two  values,  although 
a  better  agreement  can  be  obtained  by  choosing  a  higher  k  factor  for  the 
tangential  cooling.  The  largest  reasonable  value  of  k  (as  previously  dis¬ 
cussed)  gives  u*n  »  .036  Wg,  which  agrees  with  the  other_values  to  within 
experimental  uncertainties.  The  experimental  values  of  uw  tend  to  be  less 
than  theoretically  predicted  as  the  Inner  radius  Is  approached,  and  this  may  be 
a  consequence  of  the  flow  not  being  fully  developed. 

This  concludes  the  dlscussslon  of  the  mean  and  turbulent  flow  variables.  The 
flow  Is  symmetric  to  within  our  resolving  capability,  and  the  general  features 


12.  Koch,  R.  and  Feind,  K. ,  "Druckverlust  und  Warmeubergang  In  Rlngspalten," 
Chemle-Ing.-Techn. ,  30,  1958,  pp.  577-84. 
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of  the  flow  agree  with  our  knowledge  of  turbulent  shear  flow.  The  mean  velocity 
profile  Is  logarithmic,  while  the  turbulent  Intensities  have  a  parabolic  be¬ 
havior  with  minimum  values  near  the  midpoint  between  the  Inner  and  outer  radii. 
The  behavior  of  the  shear  stress  coefficient  agrees  well  with  past  experiments, 
while  uw  agrees  reasonably  well  with  the  theoretical  prediction.  Certain 
details  of  the  experiment  do  not  precisely  agree  with  theory  and  expectations 
and  this  Is  probably  a  consequence  of  the  flow  not  being  completey  developed. 


V.  TURBULENT  LENGTH  SCALES  VIA  TAYLOR'S  HYPOTHESIS 


Turbulent  length  scales  can  be  obtained  using  single  point  autocorrelation 
measurements  In  conjunction  with  Taylor's  hypothesis.  That  Is,  the 
following  correspondence  can  be  made. 


RuuUz)  -  W  Ruu(t) 
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and  W  Is  the  local  downstream  velocity.  The  correlation  functions 
are  Independent  of  z  and  t  for  fully  developed  flows.  Taylor's 
hypothesis  essentially  states  that  the  eddy  field  Is  frozen  and 
convected  downstream  by  the  local  mean  velocity.  Two  criteria 
are  necessary  for  this  assumption  to  be  valid,  (1)  the  turbulent 
Intensity  must  be  low. 


W  >>  w 

and  (2)  the  mean  flow  gradients  must  be  small  relative  to  a  typical  cross 
section  stream  eddy  size,  £ 


1  » 


£  dW 
¥  dr 


These  two  criteria  are  well  satisfied  In  the  central  annular  region 
where  mean  gradients  and  turbulent  Intensities  are  low. 


r  ■  .25  to  .75 
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These  criteria  are  not  as  well  satisfied  near  the  Inner  and  outer  radii,  but 
probably  do  not  Invalidate  the  use  of  Taylor's  hypothesis.  The  worst  case 
values  are 


w  =  .12  W 


and 


dW  „  ,  W 

lr  "  -Z  I 


The  autocorrelation  can  be  used  to  define  Integral  length  scales  by 


uu,z 


2  W 
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di 


where  Is  the  position  of  the  first  zero  In  Ruu(t);  this  convention  of 
Integrating  to  t0  Is  a  standard  one  and  Is  discussed  in  Hassan  et.  al. 

(Ref.  13).  The  Ayu  .  notation  is  one  In  which  the  first  two  subscripts 
refer  to  the  correlation  function  indices,  while  the  last  index  refers  to 
the  integration  direction  (downstream  in  this  case,  by  Taylor's  hypothesis). 

The  following  Integral  scales,  Auu  z,  AVv,z  and  A#w,z  were  obtained  as 
a  function  of  radius  and  are  presented  in  figure  (9).  The  behavior  of 
A„u  ,  and  Ayv  z  Is  expected,  where  the  length  scales  are  constant  away  from 
the* boundaries  and  go  to  zero  as  the  Inner  and  outer  radii  are  approached.  The 
two  scales  are  approximately  equal  In  the  central  region  where  Auu,z  =  .082R 
and  Avv  z  *  .090  Rq  .  The  relative  behavior  of  Ayu  z  and  Avv,z  as  the 
boundaries  are  approached  is  consistent  with  recent  theoretical  work  by  Lewellen 
and  Sandri  (Ref.  14).  Auu  z  decreases  more  rapidly  than  AyVj2  near  the 
walls  because  the  normal  velocity  (u)  Is  affected  more  readily  by  the  presence 
of  a  wall  than  the  parallel  velocity  (v).  Such  behavior  is  the  cornerstone  of 
the  Lewellen-Sandri  model. 


The  behavior  of  A^  z  is  more  complex  than  for  the  other  two  scales.  We  first 

note  that  A^  z  ~  §AUU  z  and  this  Is  a  common  feature  of  turbulent  shear 

flows.  The  length  scaiA  in  the  downstream  direction  for  the  velocity  receiving 
energy  by  direct  production  will  be  the  largest  scale  In  the  field  and  Is  typically 
2  to  8  times  the  various  other  length  scales.  The  decrease  In  near 

r  »  .43  Is  associated  with  the  zero  mean  gradient  and  a  similar  behavior  was 
observed  by  Sabot  and  Comte-Bellot  In  pipe  flow.  Results  from  their  work  are 
presented  in  Figure  10,  and  Is  *  minimum  at  r  *  0  where  dW/dr  «  0. 

The  values  for  AyU  z  are  also  shown  and  agree  with  those  observed  In  the  pre¬ 

sent  study.  The  values  of  Km,z  given  by  Sabot  and  Comte-Bellot  are  much 
greater  than  those  obtained  here  and  Hassan  et.al.;  data  from  Hassan  et.al.  Is 


13.  Hassan,  H.  A.,  Jones,  B.  G.  and  Adrian,  R.  J.,  "Measurements  and  Axisym 
metric  Model  of  Spatial  Correlations  In  Turbulent  Pipe  Flow,"  AIAA 
Report  No.  79-1562,  1979. 

14.  Lewellen,  W.  S.  and  Sandri,  G.,  "Incorporation  of  an  Anisotropic  Scale 
into  Second-Order  Closure  Modeling  of  the  Reynold  Stress  Equation," 
A.R.A.P.  Tech  Memo  No.  80-11  (Submitted  to  J.  Fluid  Mech.,  1980). 
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also  included  in  Figure  10.  (We  have  carefully  reviewed  other  works  by  Sabot 
and  Comte-Bellot,  and  they  consistently  report  values  of  greater  than 

observed  by  other  workers.)  The  increase  in  r  near  r  *  .85  in  the 
present  data  may  be  related  to  the  general  parabotic  behavior  of  A*#^  in 
pipe  flow.  Precise  agreement  between  pipe  and  annulus  flow  is  not  expected,  but 
the  behavior  of  A^  z  In  the  present  study  is  consistent  with  pipe  flow  data. 

We  compare  the  A data  with  the  von  Karman  length  scale,  Ayx  ,  defined  as 
follows  for  plane  parallel  flow, 


AVK 


dW/dy 

dZW/dy2 


where  y  is  the  cross-stream  gradient  direction.  The  motivation  for  this 
definition  Is  that  away  from  flow  boundaries,  the  scale  will  be  determined  by 
local  derivatives  of  W,  and  a  rational  choice  is  the  first  and  second  deriva 
tives.  An  analogous  expression  for  Ayx  in  cylindrical  coordinates  is 


_ dW/dr 

(l/r)(d/dr)(rdW/dr) 


Ayx  was  computed  using  our  mean  velocity  field  and  the  results  are  presented 
In  Figure  11.  The  von  Karman  length  possesses  a  number  of  features  which  are 
consistent  with  the  A^^  data,  although  these  features  tend  to  be  exaggerated. 
The  general  level  of  Ayx  Is  of  the  same  order  as  A^.z,  Ayx  is  greater  near 
the  outer  wall  than  near  the  Inner  wall,  and  Ayx  is  small  in  the  region  of 
small  mean  gradient.  It  is  remarkable  that  given  the  simplicity  of  the  von 
Karman  scale.  It  behaves  as  similarly  to  the  data  as  observed. 

The  autocorrelation  results  can  be  compared  to  concepts  derived  for  isotropic 
and  homogeneous  turbulence,  and  although  the  annulus  flow  is  neither  isotropic 
nor  homogeneous,  the  data  do  have  certain  features  In  common  with  such  flows. 

The  correlation  tensor  may  be  completely  specified  by  two  scalar  functions 
f  and  g,  for  Isotropic,  homogeneous  turbulence,  and  these  functions  are  de¬ 
pendent  only  on  the  magnitude  of  the  separation  vector.  This  Is  discussed  In 
The  Theory  of  Homogeneous  Turbulence  by  Batchelor,  (Ref. 15)  and  Figure  12  is  used 
to  define  the  f  and  g  functions.  A  separation  vector,  r,  is  specified  and 
velocities  parallel  and  normal  to  r  are  defined  as  up  ana  un.  respectively. 
The  following  definitions  are  taken? 

up(x)  up(£  +  Li  =  u2f(r) 

un(x)  un(*  +  l)  E  u2g(r) 


15.  Batchelor,  G.  K.,  The  Theory  of  Homogeneous  Turbulence,  Cambridge 
University  Press,  1967. 
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where 


by  Isotropy.  The  functions  f  and  9  depend  only  on  the  magnitude  of  £ 
and  for  Incompressible  flows,  the  following  constraint  exists. 


9  *  f  +  (r/2)(3f/3r) 


The  behavior  of  f  and  g  Is  such  that  f(0)  *  g(0)  ■  1 ,  f  will  have  a 
long  positive  tall  for  large  r,  and  g  will  go  through  zero  or  become 
negative  for  relatively  small  r. 

The  data  Is  now  placed  in  the  context  of  f  and  g  functions.  The  sepa¬ 
ration  direction  is  A z,  so  that  w  plays  the  role  of  the  parallel  velocity, 
while  both  u  and  v  play  the  role  of  the  perpendicular  velocity.  Auto¬ 
correlation  data  from  r  *  .75  Is  presented  in  Figure  13  for  the  three 
velocity  components  and  the  autocorrelation  lag  time,  t,  is  related  to  Az 
by  Wt  *  Az.  Both  Ruu  and  RVv  behave  g-llke,  becoming  zero  at  Az/Rq  *  .3, 
while  Rww  Is  f-llke,  remaining  positive  for  separations  greater  than  1.2R0. 

A  more  sensitive  test  to  Isotropic  concepts  Involves  the  length  scales  Af 
and  Ag,  defined  by 


and  the  constraint,  g  «  f  +  (r/2)(3f/3r),  yields 

VAf  ■  °'5 

The  following  correspondences  are  made, 

Ag  Auu,z  *  ^vv.z 

Af  *  Vz 

and  the  ratios  Agu,z/A|«  »  wd  Aw,z/Aww,z  are  tabulated  below  for  three 
radial  positions,  r  ■  .2/,  .43,  and  .75  where  all  lengths  are  nomallzed  by  Rq. 
(The  tllda  Is  placed  over  A  to  emphasize  that  these  scales  were  obtained  by 
Integrating  the  correlation  functions  out  beyond  to). 
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Auu,z^Aww,z 

Avv,z^Aww,z 

r  «  .27 

CVJ 

CM 

.26 

r  •  .43 

.26 

.30 

r  *  .75 

O 

CVI 

• 

CVJ 

CM 

• 

Ml  ratlos^are  well  below  the  0.5  value  associated  with  Isotropic  flow.  The 
ratio  of  Ayv.z/jLw.z  at  r  *  .43  Is  the  closest  to  0.5  because  the  flow  Is 
nearly  Isotropic  locally,  no  local  production  of  w  occurs,  and  the  v  velocity 
Is  not  as  sensitive  to  boundaries  as  u.  It  Is  not  surprising  that  the  above 
ratios  differ  substantially  from  0.5  because  of  the  use  of  Ki»%z-  The  w 
velocity  receives  energy  by  direct  production,  and  quantities  related  to  w  are 
not  expected  to  satisfy  Isotropic  relations.  In  the  following  section, 
two-point  measurements  give  a  separation  In  the  radial  direction  and  the  corre¬ 
lation  functions  for  u  and  v  are  f  and  g-like,  respectively.  Results 
for  that  case  will  be  in  closer  agreement  with  isotropic  concepts  and  the  data 
In  the  present  section  Is  given  for  future  comparison  with  the  two-point  data. 


VI.  TURBULENT  LENGTH  SCALES  VIA 
TWO-POINT  CORRELATIONS 


Two-point  correlations  for  u,  v  and  w,  and  two-point  cross-correlations 
for  u  and  w  were  measured  for  the  case  of  pure  radial  separations. 

Three  fixed  probe  positions  were  chosen,  rp  *  .27,  .43,  and  .75,  where  the 
mean  velocity  gradient  was  positive,  zero,  and  negative,  respectively.  Such 
measurements  allowed  an  examination  of  the  effects  of  mean  shear  and  boun¬ 
daries  on  the  correlation  functions.  Data  from  the  two-point  correlation 
measurements  are  presented  In  Figures  14,  15,  and  16  for  the  three  fixed 
probe  positions;  the  notation  Is  one  In  which  curves  labeled  u  refer  to 
Ruu.  defined  as 


u2  Ruu(rF;  Ar) 


The  curves  labeled  v 
similarly. 


=  11m 

J-»  OD 


/*+T 

/  u(rp,t)  u 
-T 


(rp  +  Ar,t)  dt 


and  w  refer  to  Rvv  and  R^  which 


are  defined 


Data  from  rp  «  .27,  where  dW/dr  Is  positive,  are  presented  In  Figure  14. 

In  the  region  between  the  Inner  tube  (Rj  *  .118)  and  rp,  R^  >  Ru„  >  Rvv, 
and  Rvv  Is  negative  for  r  <  .18.  The  region  outward  from  rp  snows  a 
similar  behavior,  R^  >  RUu  >  Rvv  out  to  r  »  .38;  beyond  that  point, 

Ruu  and  Rvv  possess  long  positive  tails  and  Rwm  goes  negative.  Data 
from  the  zero  mean  gradient  position,  rp  *  .43,  are  shown  In  Figure  15,  and 
the  following  relationships  hold  for  all  separations,  RUu  >  Rww  >  Rvv. 

This  result  differs  from  the  previous  case  (and,  as  will  be  seen,  from  the 
rp  =  .75  case)  and  Is  related  to  the  zero  mean  gradient.  Figure  16  shows 
data  from  rc  *  .75  where  dW/dr  was  negative.  In  the  region  approaching 
the  outer  wall  (.75  <  r  <  1.0),  R^  >  Ruu  >  Ryv»  and  Rvv  goes  negative 
for  r  i  .87.  This  behavior  Is  similar  to  that  observed  In  the  rF  *  .27 
case  for  the  region  between  the  Inner  tube,  and  rp.  This  similarity  results 
from  the  fact  that  both  regions  possess  relatively  constant  shear  and  are 
near  boundaries.  The  results  for  r£  .75  show  Rim  >  Ruu  until  r  =  .55, 
Ruu  »  Rww  for  r  s  .55.  RUu  possesses  a  long  positive  tall,  and  this  Is 
similar  to  the  behavior  In  the  rp  ■  .27  case.  In  both  instances,  RUu 
remains  positive  and  decays  slowly  as  one  moves  away  from  the  near  boundary. 
Both  Rmm  and  Rvv  go  negative  for  r  <  .75,  but  Ryy  becomes  negative  for 
much  smaller  separations  than  does  R*w. 


Cross-correlations  for  u  and  w  were  measured  and  defined  as 

r+T 

(^)V2(^)l/2  Rw(r;Ar)  =  11m  jy  /  u(r,t)  w(r  ♦  Ar,t)  dt 

and  »+T 

(^)V2(^)l/2  u^(P;AP)  =  X.  I  w(r,t)  u(r  ♦  Ar,t)  dt 

-T 


We  note  that  for  Ar  *  0, 

Ruw  (r;0)  «  R^y  (r;0)  ■  Shear  Stress  Coefficient 

while  for  Ar  f  0,  Ruw  and  Rwu  need  not  be  equal.  Data  from  the 
three  fixed  probe  positions,  rp  -  .27,  .43  and  .75  are  presented  In 
Figure  (17),  (18)  and  (19),  respectively.  The  cross  correlations  at 
rF  ■  .27  behave  such  that  Ruw  =  R*u  the  region  between  the  inner  tube 
and  rp,  where  the  mean  velocity  gradient  Is  fairly  constant,  while  for 
r  >  rF,  Rwu  >  Ruw  In  a  region  where  the  mean  gradient  Is  rapidly  decreas¬ 
ing.  Figure  (18)  shows  data  taken  about  the  zero  gradient  position,  r  *  .43, 
and  we  find  that  |RUw|  >  | Rwu I  for  all  radial  separations.  We  note  that 
Ruw(Ar  »  0)  ■  Rwu(Ar)  f  0  because  the  fixed  probe  was  not  precisely  at  the 
zero  mean  gradient  position.  It  Is  most  likely  that  RMU(Ar)  would  have  been 
zero  for  all  separations  If  the  fixed  probe  had  been  positioned  at  the  pre¬ 
cise  zero  stress  point.  The  sign  of  Ruw  Is  consistent  with  the  sign  of  the 
turbulent  stress  Ruw(r;Ar  *  0)  and  |Ruw|  reaches  a  maximum  value  of 
approximately  .18  for  separations,  Ar,  of  about  .14.  The  rp  *  .75  case 
Is  presented  In  Figure  (19),  and  the  results  are  very  similar  to  those  at 
rc  *  .27.  Rgw  -  Rwu  the  region  approaching  the  other  wall  (.75  <  r  <  1.0) 
where  the  mean  velocity  gradient  Is  approximately  constant.  In  the  region 
where  the  mean  gradient  Is  rapidly  changing,  r  <  .75,  Rwu*  >  Ruw*  and  this  Is 
similar  to  the  behavior  observed  at  rp  «  .27. 

The  correlation  functions  described  above  were  Integrated  along  Ar  to 
obtain  turbulent  length  scales  In  the  following  manner, 

/Ar 

R  (rr;Ar)d(Ar) 
uu  F  *  '  '  ' 

o 

The  notation  Is  the  same  as  discussed  In  Section  5  and  Ar0  Is  the  position 
at  which  Rab  becomes  zero.  The  various  length  scales  are  tabulated  In 

Table  I,  where  the  results  from  Section  5  are  also  Included.  Separations 
of  plus  and  minus  Ar  were  measured  and  Integral  scales  for  both  cases  are 
tabulated  and  denoted  by  the  +  and  -  signs.  All  lengths  are  normalized 
by  R0,  and  the  negative  values  for  Ayw,r  and  A*u,r  reflect  the  negative 
value  for  the  corresponding  correlation  functions. 

The  tabulated  length  scales  can  be  compared  to  pipe  flow  results  obtained 
by  Hassan  et.  al.,  Sabot,  Renault  and  Comte-Bellot  (Ref.  16),  and  Sabot  and 
Comte-Bellot.  In  all  three  papers,  a  complete  map  of  Ruu*  Rww*  Ruw*  and 
Rwu  was  obtained  for  arbitrary  values  of  Ar  and  Az  ,  and  such  data  allow 
tne  computation  of  Isocorrelatlon  curves  for  the  various  tensor  components. 
Hassan  et.  al.  tabulate  various  length  scales  at  rp  *  .65Rq  (or  .35RC  from 
the  pipe  wall),  while  Sabot  et.al.  provide  Isocorrelatlon  curves  at  r*  .50Rq  , 
but  do  not  tabulate  the  various  length  scales.  Data  from  the  present  work, 
Hassan  et.  al.  and  Sabot  et.  al.  are  tabulated  In  Table  II. 

16.  Sabot,  J.,  Renault,  J.  and  Comte-Bellot,  6.,  "Space-Time  Correlations 
of  the  Transverse  Velocity  Fluctuation  In  Pipe  Flow,"  Phys.  Fluids  16, 

No.  9,  1973,  pp.  1403-05. 
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Table  II 


Aww,r 

Aiu,r 

Vr 

Avw.r 

A*w,z 

Auu,z 

Present 
r  «  .75R0 

.100 

.093 

.051 

.079 

.242 

.082 

Hassan  et.al. 
r  ■  .65R0 

.134 

.083 

.057 

.051 

.587 

.083 

Sabot  et.al . 

*  .068 

^  .110 

1.000 

.130 

r  «  .50Rrt 
o 


All  separations 
In  direction 
away  from  pipe 
wall 


All  lengths  are  divided  by  Rp  and  are  derived  from  correlations  with  a 
radial  separation  away  from  tne  wall.  Sabot  et.al.  gave  values  for 
A^z  and  AUUf2.  while  AuM,r  and  A*u>r  were  obtained  by  reading  the 
respective  Isocorrelatlon  curves,  plotting  the  data. and  Integrating  the 
correlation  curve.  This  Is  not  a  very  accurate  technique  and  we  note  this 
by  use  of  the  approximation  symbols.  A  comparison  between  the  present 
results  and  Hassan's  shows  that  AuU  r,  AUu,z.  and  r  all  agree  to  with¬ 
in  approximately  10%,  while  A^  r  differs  by  about  30%.  Inspection  of 
Sabot's  results  shows  that  all  tfielr  length  scales  are  greater  than  observed 
here,  or  by  Hassan  and  this  was  discussed  In  Section  5.  The  difference 
In  observed  by  Sabot  et.al.  and  Hassan  et.al.  Is  disturbing,  as 

both  measurements  were  made  at  downstream  distances  of  90  diameters,  which 
Insures  well  developed  flow.  The  relatively  low  value  obtained  here  may 
be  a  consequence  of  the  different  geometry  (annulus  versus  pipe)  and  the 
lack  of  a  fully  developed  flow  field.  We  can  look  to  the  work  of  Laufer 
for  a  third  value  of  Amw,z  to  compare  with  Sabot's  and  Hassan's  values. 
Laufer  does  not  give  A iiw.z  values  directly  and  one  must  interpret  the 
given  power  spectra  and  estimate  A#*,*;  this  was  done  and  the  values  ob¬ 
tained  from  Laufer  support  the  A#w,z  ^  0.6Ro  value  of  Hassan. 

Substantial  differences  exist  for  the  behavior  of  the  cross-correlation 
functions,  Ruw  and  Ruu,  for  the  three  sets  of  data  (Sabot,  Hassan,  and 
present).  Both  Sabot  et.al.  and  we  found  | R«vu]  >  lRuwl  f°r  Pure  radial 
separations  away  from  the  outer  wall,  and  we  will  subsequently  show  that 
such  behavior  Is  consistent  with  the  production  terms  in  the  cross-correlation 
equations.  Hassan  et.al.  found  that  |Ryw|  >  |R*U|  for  pure  radial  sepa¬ 
rations,  and  the  difference  between  Sabot  et.al.  and  Hassan  et.al.  Is  again 
perplexing,  given  the  similarity  of  their  experiments. 

Me  end  this  discussion  by  noting  that  a  number  of  differences  exist  in  the 
literature  regarding  turbulent  pipe  flow.  In  Section  4,  the  question  of 
Isotropy  at  pipe  center  was  raised  and  a  careful  reading  of  Laufer,  Hassan 
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et.al.,  and  Sabot  et.al.  shows  comparable  discrepancies  between  the  ratio 
u/w  at  similar  radial  positions.  Differences  In  the  various  turbulent 
length  scales  exist  (most  notably  In  Am,,z)  and  differences  In  the  struc- 
ture  of  the  correlation  functions  (the  behavior  of  RUM  and  RWu)  are 
also  observed.  This  state  of  affairs  complicates  the  work  of  turbulent 
modelers.  The  correct  modeling  of  the  "tendency  to  Isotropy"  term  In 
2nd  order  closure  schemes  Is  very  Important  and  requires  good,  consistent 
data  on  the  relative  behavior  of  u,  v,  and  w.  As  turbulent  modeling 
becomes  more  sophisticated  and  length  scale  structure  Is  Introduced,  high 
quality,  consistent  data  on  the  various  turbulent  scales  will  be  required. 

It  Is  for  these  reasons  that  the  discrepancies  between  the  various  pipe 
flow  experiments  are  very  disturbing. 

The  two-point,  radial  separation  data  can  be  Interpreted  using  Isotropic  and 
homogeneous  turbulence  concepts  In  a  manner  similar  to  that  discussed  In 
Section  5.  The  separation  direction  was  radial  so  that  u  plays  the  role 
of  the  parallel  velocity,  while  v  and  w  play  the  role  of  the  transverse,  or 
normal,  velocity.  It  will  be  shown  that  the  w  correlations  are  g-llke 
only  at  the  zero  mean  gradient  position  where  no  direct  production  of  w 
Is  occurlng.  The  data  from  rp  ■  .27  show  that  Ruu  >  RVy  In  accord  with 
the  f  and  g-llke  character  of  Isotropic  turbulence.  For  Rj  <  r  <  .27 
Rvv  has  a  negative  region,  while  RVy  remains  positive  for  r  >  .27.  The 
behavior  of  A**  Is  not  consistent  with  Isotropic  concepts  as  R^  >  Ruu 
for  Ro<r<  .48,  and  this  Is  a  consequence  of  the  direct  production  of  w 
by  the  mean  gradient.  The  correlation  functions  at  the  zero  gradient  posi¬ 
tion,  rp  ■  .43,  behave  In  a  manner  expected  from  isotropic  theory.  The  Ruu 
correlation  Is  f-llke  and  Is  always  greater  than  Rvv  and  Rww>  while  both 
Ryv  and  Rmm  have  negative  regions  In  a  g-llke  manner.  At  rp  *  .75, 

Ruu  >  Rvv  the  two  correlation  functions  behave  f  and  g-llke,  respec¬ 
tively.  Again,  Rww  Is  not  g-llke  because  the  non-zero  mean  gradient 
results  In  direct  production  of  w. 

The  f  and  g  length  scales  can  be  defined  by  the  correspondence. 


Af  Auu,r 
Ag  Avv,r  • 

where  the  tllda  again  Is  used  to  emphasize  that  these  length  scales^are 
obtained  by  Integrating  beyond  Arp.  The  following  ratios,  Ayy,r/Auu,r 
and  Xww  r/^uu ,r»  are  tabulated  below  for  the  three  fixed  probe  positions, 
end  we  distinguish  between  positive  and  negative  Ar  by  +  and  -. 
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Avv 

,r/Auu,r 

uu 

♦ 

.72 

.70 

27 

- 

.67 

.79 

♦ 

.53 

.69 

43 

- 

.70 

.70 

+ 

.58 

1.30 

75 

- 

.47 

.90 

The  values  here  are  much  closer  to  the  ^0.5  value  than  observed  In  Section 
5,  and  this  Is  a  consequence  of  using  Auu,r»  which  Is  derived  from  a  velo¬ 
city  component  not  receiving  energy  by  direct  production.  The  values 
closest  to  0.5,  occur  at  rp  »  .43  where  the  mean  gradient  Is  zero,  and 
for  positive  separations,  which  are  In  the  direction  of  the  far  wall  and 
affected  less  by  boundaries  than  the  negative  separation  case. 

The  equations  for  the  various  correlation  functions  may  be  readily  derived 
for  the  present  flow  configuration,  and  the  derivation  Is  presented  In 
Appendix  B.  All  the  equations  have  a  similar  form,  and  we  present  the 
Rww(lL;Ajr)  equation  as  a  model. 


■ 


where  the  subscripts  1  and  2  refer  the  positions  r  and  r  ♦  Ar , 
respectively.  The  tllda  denotes  the  correlation  functions  have  not  been 
normalized  by  the  rms  velocities  (see  Section  2).  The  A mm  term  re¬ 
presents  the  advectlon  of  and  Is  equal  to 

*ww  "  NzJ/uu  +  (^)2^uu  6 

The  bracketed  term  In  Eq.  (6-1)  Is  the  production  term,  and  Is  the  one  on 
which  we  will  concentrate.  The  Dm  term  represents  turbulent  diffusion 
and  Is  comprised  of  third  order  velocity  correlations,  while  the 
term  Is  the  velocity-pressure  gradient  correlation  and  V^w  represents 
the  viscous  diffusion  term.  The  equations  fbr  the  other  non-zero  tensor 
components  are  presented  below  where  the  various  terms  are  analogous  to 
those  described  above. 


64 


The  equations  for  RUy  and  Rvv  do  not  have  direct  production  terms, 
and  the  differences  observed  for  these  functions  within  the  flow  field 
are  due  to  geometrical  effects  (boundaries)  and  the  nature  of  the  turbu- 
lent  diffusion  and  velocity- pressure  gradient  terms  (Dij  and  Pij). 

Equation  (6-1)  for  R* w  shows  a  production  term  which  provides  a  qualita¬ 
tive  explanation  of  the  behavior.  For  the  case  of  pure  axial 

separation,  (3W/3r)i  ■  (3W/3r)2,  and  at  rf  •  .43,  the  mean  gradient  Is  zero; 
this  Is  consistent  with  the  diminishing  of  Aww,z  near  the  zero  gradient 
position.  The  equations  for  Ruw  and  R^  [Eqs.  (6-4)  and  (6-5)]  show 
production  terms,  Ruu(3W/3r)2  and  Ruu(3W/3r)i,  respectively.  In  regions 
of  relatively  constant  mean  gradient  where  (3W/3r)i  »  (3W/3r)2,  one  can 
expect  Ryu  «  Rwu,  and  this  was  found  to  be  true.  Ue  write  Eqs.  (6-4)  and 
(6-5)  as  follows. 


3Ruw 

-fr 


■  -  R. 


+  Other 
Terms 


♦  Other 
Terms 


and  we  Inspect  the  behavior  at  rc  ■  .43,  where  (3H/3r) 


The  pro 


being  a  relatively  small  quantity  around  rP  ■  .43.  The  production  term, 
-Ruu(3M/3r)2  t  0  »  In  Eq.  (6-7)  Is  non  zero  as  (SW/SrJg  t  0  for  separa¬ 
tions  about  rp  •  .43  and  this  agrees  with  the  observed  behavior,  |R  |  > 
|R  |.  Also,  the  observed  sign  of  R  Is  consistent  with  the  production 
term  1"  Eq.  (6-7).  Finally,  we  observe  that  the  structure  of  the  produc¬ 
tion  terms  for  R  and  A  will  drive  |R_  (  >  |R  |  in  the  regions  of 
rapidly  changing  shear*  r  >  r_  ■  .27  and  r  <  r-  •  .75,  as  was  observed. 


VII.  SUMMARY  AND  RECOMMENDATIONS 


The  present  experimental  study  has  provided  a  complete  description 
of  the  mean  and  turbulent  velocity  fields  in  annular  flow  at  a  downstream 
distance  of  30  diameters.  A  partial  map  of  the  correlation  tensor  components 
was  also  obtained  as  well  as  various  turbulent  length  scales.  A  good  deal  of 
time  and  effort  was  devoted  to  upgrading  the  wind  tunnel  to  research  stan¬ 
dards.  and  we  are  confident  that  the  flow  field  Is  steady  and  symmetric  to 
within  such  standards.  The  downstream  position  of  30  diameters  was  Insuffi¬ 
cient  to  provide  fully  developed  flow  as  evidenced  by  the  lack  of  Isotropy 
at  the  zero  mean  gradient  position.  Even  with  this  drawback,  the  velocity 
fields  behaved  as  generally  expected.  The  mean  velocity  followed  a  loga¬ 
rithmic  profile,  while  the  downstream  fluctuations  (w)  behaved  parabollcally 
as  a  function  of  radius  with  the  minimum  occurring  at  the  zero  mean  gradient 
position.  The  relative  magnitudes  of  u  and  v  were  tn  line  with  previous  shear 
flow  results  In  the  non-zero  shear  regions,  where  w  >  v  >  u.  The  behavior  of 
the  cross  correlation  coefficient  was,  again.  In  accord  with  previous  pipe 
flow  results,  where  uw  /  |u|jw|  ^  +,  .4  In  the  regions  of  positive  and  neg¬ 
ative  mean  gradients.  The  turbulent  Reynolds  stress  was  In  approximate  agree¬ 
ment  with  the  theoretically  derived  expression  and  discrepancies  between 
theory  and  experiment  were  probably  due  to  a  lack  of  fully  developed  flow. 

The  various  turbulent  length  scales  measured  for  the  annulus  flow  were  for 
the  most  part  In  accord  with  the  results  of  other  workers.  Typically  length 
scales  for  velocity  correlations  with  pure  radial  separations  were  all  of 
order  one  tenth  the  outer  pipe  radius.  The  behavior  of  AUu,z  and  Avv.z  as 
a  function  of  radius  agreed  with  pipe  flow  results  and  their  relative  behavior 
near  flow  boundaries  provided  support  for  the  two-scale  turbulent  models  of 
Lewellen  and  Sandrl.  The  Aww.z  scale  was  a  minimum  at  the  zero  mean  gradient 
position.  The  magnitude  of  Amt,z  was  smaller  than  generally  reported  for 
pipe  flow  and  this  may  be  a  consequence  of  the  different  geometry  (annulus 
versus  pipe)  or  the  lack  of  fully  developed  flow.  The  length  scale  results 
were  interpreted  using  a  number  of  theoretical  concepts.  One  was  the  von 
Karman  length  scale  which  is  based  on  the  ratio  of  local  mean  velocity  deriva¬ 
tives.  This  simple  model  gave  surprisingly  good  agreement  with  the  general  be¬ 
havior  of  Amw ,z*  The  theoretical  framework  derived  for  Isotropic  turbulence 
was  used  to  Interpret  certain  length  scale  results.  Me  found  reasonable  agree¬ 
ment  between  theory  and  experiment  for  length  scales  associated  with  the  u  and 
v  velocities  which  do  not  receive  energy  by  direct  production.  The  best  agreement 
occured  for  length  scales  measured  about  the  zero  mean  gradient  position  where 
even  w  receives  no  energy  by  direct  production.  Finally,  the  full  equations 
for  the  correlation  functions  were  presented,  and  the  production  terms  were 
examined;  these  terms  were  found  to  be  consistent  with  our  experimental  results. 


The  present  work  has  provided  data  on  turbulent  structure  In  annular  flow. 

As  previously  stated,  much  effort  was  spent  upgrading  the  facility  to  research 
standards,  and  we  now  have  a  research  quality  apparatus  along  with  all  the 
necessary  data  taking  hardware  and  software.  A  number  of  research  projects 
could  now  be  undertaken  for  varying  degrees  of  expense  and  we  will  outline 
them.  The  first  project  would  be  an  exact  repeat  of  the  measurements  de¬ 
scribed  here  at  a  downstream  distance  of  45  diameters  (and  possibly  15  di¬ 
ameters).  This  would  provide  data  on  the  evolution  of  various  flow  quantities 
as  a  function  of  downstream  distance.  Such  Information  would  be  of  Interest 
because  not  all  flow  variables  come  Into  equilibrium  simultaneously.  Certain 
variables  obtain  equilibrium  markedly  faster  than  others,  such  as  the  Reynolds 
stress  coefficient.  This  was  evidenced  In  our  present  results  where  the  stress 
coefficient  agreed  with  well  developed  pipe  flow  results,  while  u  and  v 
were  still  evolving  (the  lack  of  Isotropy  at  the  zero  gradient  position).  Such 
a  project  would  require  no  capital  expenditures  and  could  be  done  In  approxi¬ 
mately  %  man-year.  A  study  which  could  be  done  with  a  small  capital  expenditure 
Is  a  map  of  the  various  Isocorrelatlons  at  positive,  negative  and  zero  mean 
gradient  positions.  This  would  require  correlation  measurements  for  simultaneous 
non-zero  radial  and  axial  separations,  and  a  modification  of  probes  and  trans- 
verslng  mechanisms  would  be  necessary.  Such  a  study  would  provide  a  very  com¬ 
plete  determination  of  the  effect  of  shear  on  turbulent  structure,  and  would 
again  require  approximately  h  man-year. 

The  Introduction  of  swirl  to  this  flow  field  Is  still  a  goal  of  our  re¬ 
search.  Very  little  data  exists  on  the  behavior  of  swirling  flows  and  given 
the  large  number  of  situations  where  swirl  occurs,  the  Importance  of  under¬ 
standing  such  flows  goes  without  saying.  Two  possibilities  exist  for  the  In¬ 
troduction  of  swirl  In  the  present  apparatus.  One  Is  to  run  the  tunnel  in  Its 
present  configuration  with  the  honeycomb  and  Introduce  swirl  by  use  of  vanes  or 
"twisted  tape".  The  efficiency  of  such  a  method  would  probably  be  low,  so  that 
high  swirls  could  not  be  attained.  The  precise  design  of  such  a  system  would 
require  trial  and  error,  but  should  not  be  very  expensive.  Once  swirl  Is  In¬ 
troduced,  one  can  measure  the  swirl  decay  rate  as  a  function  of  downstream 
position,  and  this  Information  would  be  extremely  useful  for  comparison  to  the 
A.R.A.P.  2nd  order  turbulence  model.  Also,  one  can  Investigate  the  effect  of 
swirl  on  the  turbulent  structure.  Such  effects  would  be  analogous  to  stably 
stratified  flow  in  the  sense  that  both  rotation  and  stratification  can  be 
stabilizing  Influences.  The  other  way  In  which  to  Introduce  swirl  Is  to  go 
back  to  the  original  design  concepts  of  this  apparatus.  The  vanes  situated 
at  the  Inlet  Introduce  swirl  which  Is  "spun-up"  during  the  contraction  process. 

It  Is  precisely  In  the  contraction  section  where  we  believe  the  asymmetries 
are  Introduced  and  modification  of  the  contraction  section  Is  In  order.  The 
contraction  section  would  be  redesigned  to  provide  a  very  smooth  and  uniform 
contraction  for  the  purpose  of  eliminating  any  asymmetries.  The  advantage  of 
such  a  technique  Is  the  generation  of  very  high  swirl  levels,  while  the  draw¬ 
backs  Involve  the  difficulty  In  creating  symmetric  flows  In  such  a  system 
and  the  fair  amount  of  construction  required.  The  rewards  would  be  many  be¬ 
cause  detailed  characteristics  of  highly  swirled  flow  could  be  studied  and 
this  would  represent  a  true  advance  in  our  understanding  of  such  flows. 
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APPENDIX  A 


The  expression  for  the  effective  cooling  velocity,  Uc,  of  a  hot  Mire  or  hot 
film  sensor  oriented  at  an  angle  a  to  the  Instantaneous  velocity,  U,  Is  given 
below  as 


Uc  *  U[cos2o  +  k2  sin2  a] 1/2 


(A.l) 


where  k  Is  an  experimentally  determined  quantity  and  the  orientation  of  the 
flow  and  sensor  are  shown  In  Figure  A-l(a).  The  output  voltage,  e,  of  an 
anemometer- 1 1nearl zer  system  will  be  proportional  to  U,  and  we  write 


e  -  vUc 


(A. 2) 


where  y  Is  the  calibration  coefficient.  The  two  sensors  of  the  X-probe  are 
labeled  A  and  6,  and  are  oriented  at  plus  and  minus  45°  to  the  mean  flow 
direction,  U,  respectively.  The  output  voltage  for  each  sensor  from  an 
Instantaneous  velocity,  U,  oriented  at  an  angle  +  to  the  mean  flow  direction 
Is  given  below  as 


eA  -  y[(U  ♦  u')2  +  (v‘)2]l/i  [cos2(eA  +  ♦)  +  k2  sin2(eft  +  ♦)] 


[■_  li/2  r  -|l/2 

(U  +  u')2  +  (v')2J  |cos2{eB  +  f)  +  k2  sin2(hg  +  *)J 


The  primed  quantities  are  the  fluctuating  velocities  and  the  orientation  of 
sensor  and  velocity  Is  shown  In  Figure  A-l(b).  The  trigonometric  functions 
are  expanded  and  we  make  use  of  the  fact  that 


tan  + 


_vl_ 
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tan  +  «  + 
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this  leads  to  the  following  expressions. 


CA  -  Jf 
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1/2 
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(i  +  4)2  +  ftyl 
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(A. 3) 


(A. 4) 


The  above  expressions  may  be  further  simplified  by  neglecting  _second_  order 
quantities  and  expanding  the  square  root  expressions  (using  u'/U,  v'/U  «  1, 
a  «  V/ff). 
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(A.6) 


The  sum  of  eA  and  eg  is  taken  as 


S  = 


2LliB  .!»(,♦ 
2  ✓  ? 


(■*J) 


(A. 7) 


and  the  average  of  es  is 


eS  *  J|(1  +  k2) 


v'jj 

a 


(A.8) 


where  is  the  new  calibration  constant.  The  fluctuating  voltage  may  be 
written  as 


eS 


eS  '  eS 


T  U 


I 


and  the  difference  of  eg  and  e^  is  taken  as 
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(A.9) 


Second  order  terms  are  again  neglected  and  we  make  use  of  the  fact  that 
eQ  -  0;  this  yields 


(A. 10) 


The  apparent  velocity  fluctuations  are  defined  in  relation  to  the  true 
fluctuations  as  follows. 


„•  -  ii  . 

apparent  Y • 

utrue 

(A. 11) 

apparent  Y ■ 

(A. 12) 

The  above  expression  shows  that  the  apparent  transverse  fluctuations  must  be 
corrected  by  the  factor  (l+k2/l-k2)  to  yield  the  true  value. 

The  k  factor  is  determined  experimentally  by  measuring  the  flow  angle  a  and 
using  Eq.(A.l).  The  measurement  is  a  difficult  one  and  a  good  deal  of  scatter 
exists  for  the  k  values  determined  here  and  by  others  (Jorgensen,  1971).  The 
value  of  k  depends  on  a,  with  the  largest  values  of  k  occurring  at  u  *  0.  At 
o  *  45°,  we  find  k  *  .32±.05  for  the  hot  film  probes  used  in  the  present 
experiment,  while  Jorgensen  (1971)  found  k  *  .36*. 04  for  comparable  hot  film 
sensors.  Friehe  and  Schwarz  (1968)  used  a  somewhat  different  expression  for 
the  effective  cooling  velocity  than  given  in  Eq.(A.l),  and  the  correction 
factor  they  give  is  consistent  with  a  value  of  k  ■  .28  for  hot  film  sensors. 
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APPENDIX  B 


The  equations  for  the  correlation  tensor  components  can  be  derived  for  each 
component  separately  as  follows,  where  the  £uw  component  Is  chosen  as  an 
example.  The  following  relationships  are  used: 


Ruw(HiAH)  “  u(r)  w(r  ♦  Ar)  u^2 


where 


and 


r  -*■  1 
r^  +  Ar  ■*  2 


The  preceeding  definitions  are  used  In  conjunction  with  the  Navler-Stokes 
equations  for  the  turbulent  fluctuations  In  cylindrical  coordinates. 


The  above  equations  are  averaged  and  added  to  yield 

{%)t  "  ♦  »'» 


where  the  following  definitions  are  used. 


[H1  ^7  +  l,2^]Ruw 


The  term  A w  represents  the  advection  of  the  correlation  by  the  mean  velocity 
while  Dy*  represents  turbulent  diffusion  and  is  comprised  of  third  order 
correlations.  PUw  represents  the  velocity  pressure  gradient  correlations 


and  Vuw  Is  the  viscous  diffusion  term.  The  equations  for  the  other  non-zero 
tensor  components  can  be  derived  1n_a  similar  manner  and  are  given  in  Section 
6.  He  do  not  bother  defining  the  A,  D,  P,  and  V  terms  for  all  the  corre¬ 
lation  equations  because  they  are  never  used  and  the  example  presented  here 


shows  the  form  of  these  terms.  The  production  term  Is  of  Interest  and  Is 


presented  In  Section  6  for  all  the  non-zero  correlation  components. 
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Figure  2  -  Schematic  of  contraction  section  where  the  placement  of  the 
polyurethane  and  honeycomb  sections  Is  shown. 
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Figure  3  -  Schematic  of  data  acquisition  system. 
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1.2 


Figure  4  -  Mean  axial  velocity  vs.  radial 


position. 
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Figure  8  -  Turbulent  shear  stres 


Figure  10  -  Turbulent  length  scales  vs.  radial  position  constructed  for  data  of 
Sabot  and  Comte-Bellot  (Ref.  8)  and  Hassan,  et.al.  (Ref.  13). 


up(x) 


up(x)up(x  +  r)  =  u2f(r) 
un(x)un(x  +  r)  *  u2  g(r) 

Up  -  u  2n  *  u2 

Figure  12  -  Schematic  of  perpendicular  and  normal  velocities  \ 

and  their  orientations  to  the  separation  vector, 
r. 
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Figure  13  -  Correlations  for  the  three  velocity  components  as 
a  function  of  axial  separation. 
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Estimated  Error 


Estimated  Error 


wu 


Figure  17  -  Cross-correlations  for  u  and  w  as  a  function  of  radial  separation; 


Estimated  Error  £ 


Figure  19  -  Cross-correlations  for  u  and  w  as  a  function  of  radial  separation;  rP  *  .75. 


